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SHOCK- TUBE  DIAGNOSTICS  UTILIZING  LASER  RAMAN  SPECTROSCOPY 


This  report  documents  measurements  of  high  temperature  gas 
properties  behind  a shock  wave  utilizing  the  laser  Raman  scattering, 

LRS,  technique.  The  data  was  obtained  at  known  post shock  conditions 

in  order  to  calibrate  the  experimental  equipment  for  gasdynamic 

laser  cavity  diagnostics  which  were  to  follow.  Since  the  LRS 

method  is  not  as  well  known  as  other  measurement  techniques,  a 

rather  complete  description  of  the  required  spectroscopic  theory  j 

is  also  presented. 
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INTRODUCTION 


The  development  o£  non-perturbing  optical  gas  diagnost;  techniques 
has  received  much  attention  in  the  past  ' • years.  This  r.  . wed 
interest  is  mainly  due  to  the  advent  o_  modern  high  powered  lasers 
and  their  associated  electro-optical  equipment.  One  such  diagnostic 
method  of  current  interest  is  the  laser  Raman  scattering  (IRS') 
technique. 

Several  experiments  have  been  reported  in  which  the  LRS  method 
was  used  to  determine  the  thermochemical  state  of  a test  gas.  The 
number  density  of  various  species  was  measured  at  low  temperatures 
in  both  static  and  flowing  gases.1  3 High  temperature  Raman  scattering 
measurements  were  also  obtained  for  temperatures  near  1700°K  in 
H2/02-air  flames.4  Several  practical  applications  of  LRS  include  the 
measurement  of  N2  Raman  scattering  in  an  arc^-heated  wind  tunnel 
environment5  and  pollutant  density  measurements  in  a pas  turbine 
exhaust. 


Hartley,  D. , "Transient  Gas  Concentration  Measurements  Utilizing 
Laser  Raman  Spectroscopy,"  AIAA  Journal  10,  pp.  687-689,  May  1972 

2 

Kellam,  J.  M.  and  Gliek,  M.  M. , "Gas  Density  Measurements  in 
a Jet  Using  Raman  Scattering,"  AIAA  Journal  ID.  no.  1389-1391. 

October  1972  1 


3 

. , ldhopf,  F.  and  Lederman,  S.,  "Specie  Concentration  Measurements 
Utilizing  Raman  Scattering  of  a Laser  Beam,"  AIAA  Journal  9, 
pp.  309-316,  February  1971 


4 

Lapp,  M. , Goldman,  L.  M.  and  Penney,  G.  M. , "Raman  Scattering  from 
Flames,"  Science  175.  pp.  1112-1115,  March  1972 

Boiarski,  A.  A.  and  Daum,  F.  L. , "An  Application  of  Laser-Raman- 
Spectroscopy  to  Thermochemical  Measurements  in  an  Arc  Heated  Wind 
Tunnel  Flow,"  Aerospace  Research  Lab.  Report  ARL  72-0126,  December  1972 

Brespwar,  G.  E.  and  Leonard,  D.  A.,  "Measurement  of  Gas  Turbine 
Exhaust  Pollutants  by  Raman  Spectroscopy,"  AIAA/SAE  9th  Propulsion 
Conference  Paper  No.  73-1276,  Las  Vegas,  Nevada,  5-7  November  1973 


NSIVC/WOL/TR  7 S - 5 3 


The  success  of  these  experiments  has  prompted  interest  in  the 
application  of  LRS  as  a diagnostic  tool  for  shock  tube  research. 

White7  suggested  using  Raman  scattering  as  a possible  means  of  making 
a noninterfering  point  measurement  of  vibrational  and/or  rotational 
temperatures  behind  an  incident  shock  wave.  He  concluded  that  the 
species  concentrations  of  individual  components  in  a gas  mixture  could 
be  measured.  Also,  these  species  measurements,  as  well  as  energy 
mode  population  measurements,  could  be  made  whether  or  not  thermo- 
chemical equilibrium  existed. 

An  attempt  was  made  by  Milling8  to  use  Raman  scattering  as  a 
diagnostic  tool  for  measuring  CO2  dissociation  rates  in  a shock  tube 
experiment.  No  results  were  obtained,  however,  due  to  the  fact  that 
the  CO2  scattering  cross  section  which  he  assumed  was  a factor  of 
approximately  ten  times  greater  than  a more  recent  measurement 
indicates . 9 

It  was  the  objective  of  the  present  study  to  observe  Raman  scattering 
from  high  temperature  02  and  N2  behind  an  incident  shock  wave  in  air. 

The  postshock  gas  temperature  and  density  were  determined  from  the 
Raman  technique  and  compared  with  known,  shock  tube  values  over  a 
relatively  wide  range  of  conditions.  Besides  this  valuable  high 
temperature  calibration,  the  experiment  also  investigated  the 
applicability  of  the  LRS  technique  to  shock  tube  excitation  time 
measurements  and  general  high  temperature  gas  diagnostics. 


7 White,  D.  R.,  "Shock  Tube  Studies  of  Nitrogen  Vibrational 
Relaxation  and  Methane  Oxidation,"  Aerospace  Research  Lab.  Report 
ARL  70-0107,  June  1970 

Milling,  R.  W. , "Use  of  Time  Resolved  Raman  Scattering  for  the 
Measurement  of  Dissociation  Rates  in  a Shock  Tube,"  Seiler  Research 
Lab.  Report  SRL  69-005,  June  1969 

g 

Penney,  C.  M. , Goldman,  L.  M.  and  Lapp,  M. , "Raman  Scattering 
Cross  Sections,"  Nature  Physical  Science  2.35,  No.  58,  pp.  110-112, 

7 February  1972 
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ri.  DETERMINATION  OF  THE  TIIERMOCIIEMICAL  STATE  OF  A SHOCK  HEATED  CAS 


A.  GENERAL  DISCUSSION  ] 

When  light  passes  through  a gas,  weak  randomly  directed  radiation 
at  the  same  frequency  will  appear  due  to  scattering  of  the  light  by 
the  gas  molecules.  The  blue  of  the  sky  represents  such  an  effect. 

In  1928  Raman10  observed  the  spectrum  of  the  light  scattered  from 
many  different  gases  and  liquids.  Using  a monochromatic  source  of 
frequency,  v , he  noted  that  the  scattered  light  not  only  contained 
this  incident  frequency  (i.e.,  Rayleigh  scattering)  but  also  several 
shifted  frequencies,  v0-vj,  Vq-v^,  and  others,  and  v0+v£.  It 

was  found  that  v-^  and  V2  were  frequencies  which  were  characteristic  of 
the  scattering  molecule.  That  is,  for  a particular  molecule,  the 
shifted  frequencies  were  different  than  frequencies  obtained  for  any 
other  molecule.  This  frequency  shift  phenomenon  became  known  as  the 
Raman  effect.  The  radiation  shifted  toward  lower  frequencies  is 
called  "Stokes"  radiation  and  the  weaker  radiation  at  higher  frequencies, 
"anti-Stokes."  Since  its  discovery,  the  Raman  scattering  effect  has 
been  examined  for  many  different  molecules  and  has  been  an  important 
tool  used  by  physicists  for  studying  molecular  structure  at  known 
densities  and  temperatures. 

The  Raman  phenomena  can  be  qualitatively  explained  by  considering 
the  molecule-photon  interaction  from  the  standpoint  of  conservation 
of  energy.  The  incident  photon  at  frequency  vD  and  energy  Eo^hcvo 
interacts  with  a molecule  having  energy  Em.  This  inelastic  interaction 
produces  a scattered  photon  of  a different  energy,  ERaman=hcvRaman > 
at  a shifted  frequency  VRaman.  The  energy  differential  is  accounted 
for  through  a gain  or  loss  of  energy  by  the  molecule.  Since  the 
molecular  energy  is  quantized,  only  discrete  amounts  of  energy  may  be 
exchanged  in  the  interaction.  Hence,  the  molecule  must  have  an 
energy  E^  after  the  interaction  which  differs  from  its  initial  energy 
by  a specified  amount  AEr0+.v^j).  This  energy  differential  is  a 
known  function  of  the  bound  energy  states  of  the  particular  scattering 
molecule.  Balancing  energy  before  and  after  the  interaction,  the 
photon  energies  (i.e.,  frequencies)  observed  in  Raman  scattering  must 
then  correspond  to  the  photon  energy  of  the  incident  radiation 
shifted  by  an  amount  which  is  characteristic  of  the  scattering  molecule. 


■^Raman,  C.  V.,  Indian  Journal  Physics  2,  p.  287  , 1928 
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The  amount  of  frequency  shift  is  given  by: 

AE_ 


v ± 

o Raman 


"rot-vib 

"he 


CD 


At  first  there  seems  to  be  a myriad  of  possibilities  fer  the 
scattered  photon  frequencies  due  to  the  many  possible  molecular 
vibrational  energy  states,  denoted  by  quantum  number  v,  each  having 
many  associated  rotational  energy  levels,  denoted  by  the  quantum 
number  J.  However,  there  are  selection  rules  which  govern  molecular 
energy  transitions  that  can  occur  between  various  energy  levels. 

The  selection  rules  for  a diatomic  molecule  undergoing  Raman 
transitions  are: 


.s 


and 


Av  = 0,±1  (for  vibrational  transitions) 


AJ  = 0,±2  (for  rotational  transitions) 


(2) 


These  selection  rules  limit  the  number  of  allowed  frequency  shifts 
tremendously.  However,  in  this  study  even  fewer  transitions  will  be 
investigated  since  only  the  Av=+1  vibrational  and  AJ=0  rotational 
transitions  are  important.  These  transitions  give  rise  to  the  so- 
called  Stokes  vibration-rotation  O-branch  Raman  effect. 

Even  with  the  above  restrictions,  many  transitions  are  still  possible. 
This  becomes  evident,  for  example,  when  considering  the  Stokes 
Q-branch  v=0  to  v=l  vibrational  transition,  where  J=1  to  J=l,  J=2 
to  J»2,  etc.  rotational  transitions  can  occur  as  well  as  the  J=0 
to  J=0.  Due  to  small  differences  in  rotational  energy  level  spacing, 
the  resulting  frequencies  are  quite  close  but  distinguishable.  Also, 
v=l  to  v=2,  v=2  to  v=3,  etc.  vibrational  Q-branch  transitions  can 
occur  with  all  the  AJ=0  transitions  included  in  each  Av=l  vibrational 
transition.  The  various  vibrational  transitions  are  easily  distinguished 
due  to  molecular  anharmonieity . Hence,  1000-2000  Raman  O-branch 
transitions  may  be  considered  in  the  analysis  of  the  Stokes 
rotation-vibration  spectrum. 

B.  ROTATION -VIBRATION  LINE  INTENSITIES 

Smekal11  predicted  the  Raman  effect  in  1923  and  since  that  time 
•the  theory  of  Raman  scattering  has  become  well  developed  from  both 


11 


Smekal,  A.,  "Zur  Quantentheorie  der  Dispersion, 
Naturwiss  ll,  p.  873,  September  1923 
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the  quantum  mechanical  and  the  semi-classical  viewpoints.12*13 
Hence,  only  the  results  of  the  detailed  theoretical  investigations 
will  be  presented  in  this  report. 

The  observed  Raman  intensity  flux  in  photons/second,  resulting 
from  a particular  Stokes  rotation-vibration  transition  is  given  by: 


N. 


v,  J 


* NoXv,jVj“n 


(4) 


where  NQ  is  the  number  of  incident  photons/second,  % is  the  optically 
observed  length  of  the  laser  beam  in  centimeters,  ft  is  the  solid 
angle  of  the  collection  optics  in  steradians,  n is  the  number  of 
molecules  per  cubic  cm  in  the  gas  sample  and  Xv  J is  related  to  the 
fraction  of  molecules  in  energy  state  Ev  j.  Equation  (4)  shows  that 
collected  intensities  vary  directly  with’the  incident  photon  flux. 

This  fact  explains  why  the  advent  of  modern  high  intensity  lasers 
has  increased  interest  in  Raman  scattering  as  a practical  diagnostic 
technique.  Also,  the  Raman  intensity  is  directly  proportional  to  gas 
number  density  so  measurement's  can  be  used  to  determine  unknown 
density  conditions. 

Further  examination  of  Equation  (4)  shows  that  the  collected 
intensity  is  dlirectly  related  to  several  experimental  factors,  namely, 
collector  solid  angle  and  observed  laser  beam  length.  Since  the 
Raman  photon  flux  is  scattered  nearly  uniformly  into.  4n  steradians,  a 
larger  collector  implies  a greater  collection  efficiency.  The 
scattering  volume  (i.e.,  volume  from  which  collected  nhotons  originate) 
is  defined  by  the  laser  beam  cross  section  dimensions  and  the  field 
of  view  of  the  optics.  The  volume  varies  linearly  with  the  square  of 
beam  diameter  and  observed  beam  length.  This  is  an  important  aspect 
of  the  Raman  technique  since  it  indicates  that  the  measurement  is 
truly  a point  determination  of  gas  properties  within  the  scattering 
volume.  By  focusing  the  laser  beam  and  using  optical  stops  to  obtain 
a small  observed  length,  the  spatial  resolution  can  be  improved. 

However,  Raman  signal  intensity  is  sacrificed  when  smaller  and 
smaller  field  stops  are  employed. 

Finally,  the  relation  of  the  Raman  scattering  intensity  to  energy 
state  population  implies  that  detailed  information  can  be  obtained  on 
population  distributions  in  a gas.  For  non-Boltzmann  conditions  such 
as  would  exist  in  a laser  cavity,  a population  distribution  measurement 
is  important.  However,  when  Boltzmann  statistics  prevail,  determination 
of  the  temperatures  of  various  molecular  modes  is  critical. 


12 

Kramers,  H.  A.  and  Heisenberg,  W. , "Uber  die  Strenung  von 
Strahlung  durch  Atome,”  Z.  Physik,  31,  p.  6Bl,  January  1925 

^Placzek,  G. , ’’The  Rayleigh  and  Raman  Scattering,”  IJCLR-Trans- 
526(l),  translated  from  a publication  of  the  Akademische  Verlagsgesell - 
schaft  G.M.B.H. , Leipzig,  1934 
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The  symbol  oy  j in  Eq.  C4)  is  the  scattering  cross  section 
(i.e.,  cnr /particle-sr)  or  transition  probability  for  the  particular 
v-*-v+l,  J -+J  transition  and  is  written  as: 

=v,J  - K ("6'^v.jJ4  8JSv  t5) 

where  K is  a constant  which  can  be  theoretically  estimated  and 
experimentally  measured.  The  frequency  shift  Avy  j of  a Stokes  Raman 
Q-branch  line  which  appears  in  Eq.  (5)  is  given  by  Reference  14  as:14 

Avv,J  = Avv-aeJ(J+l)-8eJ2(J+l)2  (6) 

where : 


Avv  = u,e-2a)exe^v+1^+U)eye^3v  +6v+13/4^ 


In  Equations  (6)  and  (7)  above,  u>g,  wexe  and  a)eyp  are  vibrational 
constants  and  ote,  8e,  are  rotational  constants  of  the  scattering 
molecule.  The .values  of  the  above  parameters  used  in  the  calculations 
reported  in  this  text  were  taken  from  Reference  14.  The  symbol  gj 
in  Eq.  (5)  is  the  spin  degeneracy  factor,  which  depends  on  the  symm  -try 
properties  of  a particular  molecule  and  whether  J is  an  even  or  odd 


integer, 


the  vibrational  band-strength  factor,  is  given  by: 


sv  = v+1 


From  Eq.  (8)j  it  can  be  seen  that  a v=l  to  v=2  transition  has  twice 
the  probability  of  occurring  than  the  v=0  to  v= 1 . For  a Boltzmann 
population  distribution,  xy  J is  a function  of  the  vibrational 
temperature  Tv  and  the  rotational  temperature  Ty  of  the  scattering 
molecule  and  is  written: 


S ,e 


1.44F(J0 

- T - 


1 . 44G(v' 
T 

v 


‘rot^vib 


^Herzberg,  G.,  Spectra  of ^Diatomic  Molecules,  D.  van  Nostrand 
Company,  Inc . , New  York,  pp.  552-561 , 19  50 
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where  the  state  sums  Q . and  Q . , are  given  as 

rot  vie 


<rot 


1.44F(J) 
(2J+l)e  Tr 


(9) 


Q . 

'If  1 


vib 


1.44GO) 
2 e Tv 


the  functions  F(J)  and  G(v)  which  appear  in  Eq.  (9)  are  so-called 
term  values  and  are  written  as: 


F(J)  = BvJ(J+1)-DvJ2CJ+1)2 

and 


(10) 


G(v)  = u,eCv+l/2)-uexe(v+l/2)2+weyeCv+l/2)?1 


In  Eq.  (9)  Sj  is  the  line  strength  factor  for  Q-branch  transitions 
and  is  given  by: 

Sj  = (2J+1)  J-0 , 1 , 2 * " * (11) 


In  Eq.  (10)  Bv  and  Dy  are  written  as: 


Bv  " Be  ' “e  fv+1/z5 
ny  = De  + Be  (v+1/2) 


(12) 


Again  Be  and  De  are  rotational  constants  of  the  scattering  molecule 
(Ref.  14). 


The  above  equations  show  that  the  Raman  line  intensities  are 
complex  functions  of  number  density,  rotational  temperature,  vibrational 
temperature  and  various  molecular  constants  and  known  experimental 
factors.  Further,  Raman  frequency  shifts  are  uniquely  dependent  on 
-the  scattering  molecule.  Hence,  a measurement  of  the  frequency 
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spectrum  can  yield  an  identification  of  the  molecules  which  are 
present  in  the  test  sample.  Also,  the  concentration  of  each  constituent 
and  the  rotational  and  vibrational  temperatures  of  each  molecule 
species  can  be  determined  from  a spectral  intensity  analysis. 

C.  TEMPERATURE  DETERMINATION 

In  the  above  equations  the  Raman  spectrum  (i.e.,  N„  j vs. 

^vv,j)  was  shown  to  be  a function  of  the  rotational  and ’vibrational 
population  distributions  of  the  molecules  in  a gas.  For  Boltzmann 
population  conditions,  the  equations  imply  that  a measurement  of 
the  spectrum  would  yield  information  concerning  the  rotational  and 
vibrational  temperatures  independently.  Hence,  thermodynamic  non- 
equilibrium conditions  (i.e.,  Tv  t Tr)  could  be  investigated.  The 
above  statements  come  as  no  great  surprise  to  one  versed  in  other 
spectroscopic  techniques,  for  the  power  of  spectroscopy  has  been  its 
ability  to  determine  detailed  thermochemical  information  about 
molecules  in  a test  gas.  The'  following  describes  how  Raman  spectroscopic 
information  can  be  utilized  in  this  respect. 

From  the  above  equations,  the  temperature  dependence  of  the 
Q-branch  Raman  rotation-vibrational  lines  of  molecular  nitrogen  was 
calculated  for  illustrative  purposes.  The  results  for  the  even 
rotational  lines  at  Tv  = Tr  * 300°K  and  100(fK  are  plotted  in 
Figure  1.  A shift  in  the  peak  rotational  intensity  within  the  (0,1) 
rotation-vibration  band  at  2300  cm"1  occurs  as  the  temperature 
increases.  This  so-called  "blue  shift"  is  due  to  the  increase  in 
population  of  the  upper  rotational  states  with  increasing  rotation 
temperature.  Also,  at  100CPK  a small  set  of  frequency  shifts  appear 
at  about  2300  cm'1  resulting  from  the  (vsl-»-v=2)  or  (1,2)  Q-branch 
transitions.  These  upper  state  Raman  transitions  appear  in  the  spectrum 
due  to  the  increase  in  the  population  of  the  upper  vibrational  states 
with  increasing  vibrational  temperature. 

In  practice,  a spectrum  such  as  that  shown  in  Figure  1 would  be 
difficult  to  obtain  due  to  the  high  resolution  required  to  give 
individual  rotational  line  intensities  and  due  to  the  overlap  of 
rotational  lines  of  adjacent  Q-branches.  Hence,  resolved  spectra 
will  not  normally  be  obtained  experimentally.  To  compute  the  unresolved 
spectrum,  the  Raman  line  intensities  given  above  must  be  convoluted 
with  the  apparatus  function  A(AX)  of  the  measuring  instrument  since 
more  than  one  Raman  line  will  be  transmitted  by  the  device.  The 
apparatus  function,  with  a peak  transmission  at  wavelength  A , may 
then  be  scanned  through  the  spectrum.  The  unresolved  spectral 
intensity  is  then  given  by: 


Nv,J  ACAA) 


1:6 
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AA  = Ap  - Av,J 


V’J  ' VAvv,J 


Nv>j  and  Avvj  are  the  Raman  line  intensities  and  frequency  shifts, 
respectively^  given  by  the  detailed  equations  in  the  previous  section. 

By  varying  A through  the  region  o£  the  Raman  lines,  the  unresolved 
spectral  scan  can  be  determined  theoretically  using  a digital 
computer  program  to  perform  the  straightforward  but  tedious  computations 
The  computer  code  does  nothing  more  than  take  each  of  the  few  thousand 
Raman  line  intensities  and  sum  the  product  of  these  intensities  and 
the  apparatus  function  transmission  values  at  the  Raman  wavelength  of 
the  particular  transition.  Bach  sum  corresponds  to  a particular 
wavelength  position  of  the  peak  of  the  apparatus  function  (or 
transmission  curve  if  you  will).  The  computer  then  increments  this 
peak  wavelength  by  a specified  amount  and  repeats  the  entire  procedure 
again  and  again  until  a theoretical  wavelength  scan  is  produced.  Most 
of  the  Raman  lines  will  not  add  to  the  sum  of  a narrow  apparatus 
function  because  the  transmission  is  zero  at  their  particular  wave- 
length, also  when  A is  outside  the  range  of  the  larger  Raman 
intensities,  the  sum  will  be  negligible  even  at  the  peak  transmission 
of  the  apparatus  function.  An  example  of  such  a calculation  for 
molecular  nitrogen  at  a constant  number  density  and  temperatures  of 
300°K  ando3000°K  is  sho^n  in  Figure  2 for  a ruby  laser  source  where 
A0  = 6943A.  A narrow  5A  half-width  triangular  apparatus  function 
was  assumed.  Note  that  individual  rotational  lines  are  not  detectable 
as  in  Figure  1 since  the  apparatus  function  is  wide  enough  to  include 
many  lines  at  once.  However,  the  various  Q-branches  are  able  to  be 
Tt-solved  with  this  particular  apparatus  function.  Also,  note  the 
,ast  difference  in  the  spectra  at  the  two  temperatures  shown  in 
Figure  2.  At  300°K  the  Raman  spectrum  is  limited  to  a narrow  region 
of  the  total  spectrum  while  at  3000°K  a much  broader  wavelength 
interval  is  covered. 


Since  most  early  Raman  work  was  done  at  room  temperature  or  below, 
the  term  Raman  "line"  was  adopted,  because  the  Raman  spectrum  looked 
like  a "line"  (e.g.,  see  Figure  2).  However,  approximately  20 
actual  rotational  lines  make  up  the  O-branch  spectrum  at  room  temper- 
ature, as  shown  in  Figure  1;  and,  therefore,  the  term  "line"  is 
incorrect.  At  higher  temperatures,  where  other  O-branches  appear  due 
to  upper  level  vibrational  transitions,  the  term  "lino"  becomes 
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meaningless  as  well  as  incorrect.  Therefore,  a better  term  would  be 
Raman  "rotation-vibration  spectrum"  or  just  "Raman  spectrum,"  which 
has  been  adopted  in  this  report. 


The  apparatus  function  used  in  the  calculations  plotted  in 
Figure  2 was  somewhat  narrow  for  the  spectral  analyzer  used  in  this 
study.  Therefore,  the  actual  apparatus  function  employed  in  the 
present  experiment  was  used  to  compute  the  unresolved  Raman 
spectrum.  The  results  of  these  calculations  are  shown  in  Figure  3 for 
oxygen  ai  temperatures  of  Tv  = Tr  = 300,  000  and  1500°K.  Note  that 
the  various  Q-branches  are  not  resolved  for  this  relatively  broad 
apparatus  function,  yet  the  spectral  intensity  profile  is  quite 
sensitive  t'1  gas  temperature. 


For  a particular  Av=l  transition  (i.e.,  (0,1)  O-branch), 
Eq.  (4)  can  be  written: 


0,J 

‘0,0 


VAv0,J 

vo*Av0,0 


Sjgj 

S080 


1.44FCJ) 

T 


(34) 


i 


Thus,  the  relative  line  intensity  distribution  in  a given  O- 
branch  is  a function  only,  of  the  rotational  temperature  of  the  gas. 

This  can  be  readily  observed  in  Figures  1 and  2,  inhere  the  wavelength 
of  peak  intensity  in  the  (0,1)  O-branch  has  shifted  toward  the  blue 
end  of  the  spectrum  at  the  higher  temperature  condition.  Also,  note 
that  the  asymmetry  of  the  spectral  profile  is  quite  different  for  the 
elevated  temperature  case.  Thus,  by  curve  fitting  the  spectral  profile 
of  the  (0,1)  Q-branch  and/or  any  other  branch,  the  rotational  temperature 
of  the  gas  can  be  obtained  from  the  rotation-vibration  spectrum. 

Since  the  branches  must  be  resolved,  the  narrower  apparatus  function, 
such  as  that  used  in  Figure  2,  must  be  used  in  order  to  obtain 
information  on  the  rotational  state  of  the  scattering  molecules. 

Summing  Eq.  (4)  over  the  rotational  quantum  number  J will  result 
in  the  following  expression: 


where 


N, 


v - Av 
o v 


VAvo 


v 


1 . 4 4 G ( v ) 
T 


v 


(15) 


Av  = a)  -2w  x (v+l)+o>  y (3vi'+6v+13/4)  (16) 

V 0 0 6 0 0 
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Hence,  the  relative  integrated  intensity  of  the  branches  is  indicative 
of  the  vibrational  state  of  the  gas.  In  other  words,  the  shape  of 
each  branch  is  a function  of  the  rotational  temperature,  but  the 
relative  intensity  of  the  various  branches  is  dependent  only  on  th«' 
vibrational  temperature  of  the  scattering  molecules.  It  must  be 
pointed  out  that,  in  general,  accurate  integrated  Q-branch  intensities 
are  difficult  to  obtain  due  to  the  overlap  of  the  various  branches. 
Thus,  Eq.  (14)  cannot  be  used  in  most  cases  to  calculate  the 
vibrational  temperature  from  the  unresolved  spectrum.  Instead,  curve 
fitting  techniques  must  be  employed  to  obtain  the  vibrational 
temperature  of  the  gas.  Note  that  both  of  the  above  spectral  analysis 
methods  require  only  relative  intensity  measurements;  no  absolute 
intensities  are  needed. 

For  equilibrium  conditions  (i.e.,  Tv=Tr=T) , a curve  fit  of  the 
entire  spectral  intensity  profile  can  be  used  to  obtain  the  gas 
temperature.  Either  narrow  or  broad  apparatus  functions  as  in 
Figure  2 or  Figure  3 can  be  used.  However,  if  thermodynamic 
nonequilibrium  is  suspected,  a separate  curve  fit  must  be  made  of  one 
or  more  Q-branches  to  determine  the  rotational  temperatures 
independently  of  the  vibrational  state  of  the  gas.  Then,  the 
vibrational  temperature  can  be  determined  from  a curve  fit  of  the 
entire  rotation-vibration  spectrum.  The  rotational  temperature  of 
the  gas  can  be  used  as  data  for  this  calculation.  Note,  that 
narrow  apparatus  functions,  as  in  Figure  2,  must  be  used  for  the 
nonequilibrium  case  in  order  to  determine  independently  the 
rotational  and  vibrational  states  of  the  test  gas. 

D.  NUMBER  DENSITY  MEASUREMENT 

A cursory  examination  of  Eq.  (4)  might  imply  that  the  Raman 
spectrum  could  be  used  directly  to  evaluate  the  number  density,  n, 
of  a gas  under  varying  temperature  conditions.  This  is  not  true 
since  Raman  line  intensities  are  dependent  upon  temperature  and’ 
number  densities  simultaneously.  This  coupling  of  number  density 
and  temperature  is  common  to  all  spectroscopic  techniques  and 
arises  from  temperature  effects  on  the  population  distribution  of 
the  various  energy  modes.  The  coupling  effect  can  be  readily 
demonstrated  by  considering  the  spectral  plot  shown  in  Figure  2, 
where  constant  number  density  was  used  in  the  calculations.  Note 
that  if  the  (0,1)  Q-branch  peak  intensity  were  to  be  used  to  compute 
number  density  .at  3000°K  without  making  a correction  for  temperature 
effects,  the  calculated  values,  compared  to  room  conditions,  would  be 
in  error  by  a factor  of  about  5.  Hence,  the  gas  temperature  must  be 
known  in  order  to  obtain  the  correct  number  density  of  any  of  the 
species  in  the  test  gas.  A temperature  determination  is,  however, 
not  dependent  on  knowledge  of  the  actual  density.  For  nonequilibrium 
conditions,  both  the  vibrational  and  rotational  temperatures  must  be 
obtained  in  order  to  correct  for  temperature  dependence.  This 
temperature  effect  on  number  density  measurements  can  be  minimized 
under  certain  circumstances.  Note  that  by  using  a broader  apparatus 
function,  as  shown  in  Figure  3,  the  resulting  oxygen  Raman  spectrum 
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is  nearly  temperature  independent  at  7 7 7 8 A - Hence*  by  choosing  proper 
experimental  conditions  the  need  for  temperature  corrections  can  be 
eliminated  if  a certain  amount  of  calculable  error  is  tolerable.  In 
this  case,  the  intensity  of  the  oxygen  Raman  spectrum  at  7778A  is 
nearly  directly  proportional  to  the  number  density  of  oxygen  for  the 
particular  apparatus  function  used  in  the  calculations.  Varying 
the  shape  of  the  apparatus  function  will  change  the  optimum  wave- 
length and  also  the  error  estimates. 

To  determine  the  number  density  of  the  constituents  of  high 
temperature  air,  the  measurement  of  absolute  spectral  intensities  is 
not  required.  Those  difficult  absolute  measurements  may  be  avoided 
through  the  application  of  a calibration  point  which  is  the  known 
composition  and  number  density  of  room  air  existing  before  and/or 
after  the  high  temperature  test  condition.  Also,  the  number  density 
of  any  species  which  is  present  only  in  the  test  environment  can  be 
determined  by  ratioing  its  intensity  value  to  a room  temperature 
constituent  such  as  nitrogen.  Of  course  the  nitrogen  would  also  have 
to  be  present  in  the  test  environment.  This  is  an  attractive  feature 
of  the  Raman  scattering  diagnostic  method  over  other  spectroscopic 
techniques  which  require  either  absolute  intensity  measurements  or 
complex  calibration  procedures. 

E.  RATE  MEASUREMENTS 

1.  Vibrational  Excitation  Time 

A finite  time  is  required  for  the  vibrational  temperature  to 
equilibrate  with  temperatures  of  the  rotational  and  translational  modes 
of  motion  during  rapid  changes  in  molecular  energy.  This  nonequilibrium 
process  can  be  investigated  using  Raman  scattering.  To  monitor  both 
the  vibrational  and  rotational  temperature  would  be  straightforward 
but  cumbersome.  A narrow  apparatus  function  must  be  used,  as  discussed 
above,  and  data  must  be  simultaneously  obtained  from  at  least  three 
different  wavelength  positions  requiring  a good  deal  of  instrumentation. 

This  type  of  temperature  data  would  be  informative;  however,  often 
only  the  time  for  the  vibrational  mode  to  equilibrate  (i.e., 
vibrational -excitation  time  or  ry.-p)  is  of  primary  interest.  A 
simple  method  can  be  used  to  determine  ty-x  utilizing  Raman  spectroscopy 
by  monitoring  the  Raman  intensity  using  a single  spectrum  analyzer. 

In  this  case,  the  transmission  peak  of  the  apparatus  must  be  located 
at  the  proper  wavelength  position  in  the  spectrum.  . Raman  intensity 
calculations  were  performed  for  oxygen  using  Eq.  (13)  and  assuming 
an  abrupt  change  in  the  rotational  temperature  from  ambient  conditons 
to  a value  of  1100°K,  due  to  the  passage  of  a shock  wave.  The 
vibrational  temperature  was  then  theoretically  increased  from  300°K 
to  1100°K  (i.e.,  from  ambient  to  equilibration  with  the  rotational 
temperature)  as  occurs  behind  the  wave  during  vibration  excitation. 
Results  of  these  calculations  are  shown  in  Figure  4.  Note  that  there 
is  a significant  change  in  the  spectrum  with  just  the  rotational 
mode  excited.  Then  the  spectrum  changes  .further  as  the  population 
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of  the  upper  vibrational  state  increases.  Hence,  ty.j  is  the  time 
between  the  first  change  just  behind  the  wave  and  the  attainment  of  a 
constant  value  of  the  Raman  intensity  at  thermodynamic  equilibrium. 

Note  that  in  Figure  4 the  broad  22. 5A  apparatus  function  is  sufficient 
for  these  excitation  time  experiments. 

2.  Chemical  Reaction  Time 

As  with  the  internal  modes  of  motion,  chemical  reactions  also 
take  a finite  time  to  occur.  Usually  the  excitation  time  for  the 
internal  modes  is  much  shorter  than  reaction  times  so  coupling  can 
often  be  neglected.  Measurement  of  the  chemical  reaction  rate  using 
Raman  scattering  involves  monitoring  the  formation  or  consumption  of  a 
particular  constituent  with  time.  However,  for  a pure  gas  the  chemical 
reaction  process  alters  the  heat  content  which  results  in  a temperature 
change  during  the  reaction.  As  discussed  in  Section  D,  spectroscopic 
intensities  are  functions  of  temperature  as  itfell  as  density;  Hence, 
errors  could  be  introduced  into  the  rate  measurements  if  this 
temperature-density  coupling  is  not  taken  into  account. 

Again,  one  may  monitor  both  the  temperature  and  density.  Thus* 
Raman  intensity  data  would  have  to  be  obtained  simultaneously  for  at 
least  two  different  wavelength  positions.  The  measured  temperature 
could  then  be  used  to  obtain  corrected  density  histories  from  which 
chemical  rates  could  be  computed.  A simpler  way  which  avoids  the 
temperature-density  coupling  would  involve  diluting  the  chemical 
reactants  in  a buffer  gas  such  as  argon.  The  chemical  reaction  would 
result  in  a negligible  change  in  the  total  gas  energy  so  the 
temperature  would  remain  nearly  constant.  Now,  a single  measuring 
instrument  (i.e.,  one  broad  apparatus  function]  could  be  used  to 
directly  monitor  density  changes  with  time.  Chemical  reaction  rates 
could  then  be  calculated  without  any  error  introduced  by  the  need 
to  measure  temperature. 
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III-  SHOCK  TUBE  IN  RAMAN  SCATTERING  EXPERIMENTS 


A.  GENERAL  DISCUSSION 

The  shock  tube  is  essentially  a device  in  which  a planar  shock 
wave  is  generated  by  the  sudden  bursting  of  a diaphragm  which  separates 
a high-pressure  low  molecular  weight  gas  from  a test  gas  at  low- 
pressure.  Behind  this  wave  the  experimental  gas  is  brought  almost 
instantly  to  a known  and  controlled  high  temperature,  held  at  steady 
temperature  and  pressure  for  a few  hundred  microseconds,  and  then 
cooled  suddenly.  If  the  gas  is  studied  at  some  point  along  the  tube, 
then  it  is  possible  to  follow  the  progress  of  chemical  reactions, 
and  vibrational  excitation  processes  as  the  sample  of  heated  gas 
passes  the  measurement  point.  Since  the  density  of  the  gas  behind 
the  wave  is  usually  quite  high,  then  Raman  scattering  measurements 
are  facilitated  by  the  Increase  in  scattering  intensity.  The 
generation  of  a known  thermochemical  state  of  high  density  gas 
implies  that  the  shock  tube  is  a valuable  tool  for  calibrating 
Raman  measurement  apparatus  and  validating  Raman  theory  at  elevated 
temperature  conditions. 

B.  BRIEF  DESCRIPTION  OF  FACILITY 

To  generate  the  temperature  and  density  conditions  necessary  for 
this  study,  a 25  mm  I.D.  3.5  meter  long  shock  tube  was  used.  The 
driver  was  operated  with  helium  while  the  driven  tube  contained  either 
pure  nitrogen  or  dry  air.  Scored  aluminum  diaphragms  were  used  and 
they  were  punctured  with  a hand  operated  plunger  after  preset  charges 
of  gas  were  placed  in  the  tube.  The  plunger  method  provided  repeatable 
shock  velocities  which  were  essential  to  obtain  accurate  data  for 
comparison  with  shock  tube  theory.  A small  dump  tank  was  used 
following  the  test  section. 

Driven  pressure,  pjf  was  measured  with  a 0-400  torr  Wallace- 
Tiernan  gauge  and  driver  pressure,  pj,  was  accurately  and  repeatablv 
preset  with  a pressure  transducer  and  digital  readout  system. 

The  wave  speed  was  determined  for  each  shot  using  three  pressure 
transducers  to  start  and  stop  two  Monsanto  Model  401  timers  which 
were  accurate  to  within  one  microsecond.  The  pressure  transducers 
were  located  on  either  side  of  the  test  section  and  the  shock  speed 
at  the  point  of  the  Raman  measurements  was  computed  assuming  linear 
shock  speed  attenuation.  Error  estimates  showed  that  the  computed 
speeds  were  accurate  to  within  2 percent. 
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\vas  t0  provide  a trigger  signal  to  fire 
Another  alteration  was  the  addition  of  a Raman  test  chamber 
the  tube.  Finally  a support  was  attached  to  the  end 
...  bolted  to  the  floor.  The  purpose  of  this  stiffener  was 

vibration  during  the  diaphragm  rupturing.  Removing 

^oV^ra5*°n  eliminated  any  significant  misalignment  of  the  beam  with 
the  collection  optics. 

C . P OS  TS  H OC  K P R OP  E RT I ES 

A pressure  limit  on  P4  of  325  psia  was  arrived  at  from  safety 
considerations  since  the  diaphragm  plunger  was  hand  operated.  Also 
to  obtain  reasonably  high  postshock  densities  pi  was  kept  above 
b0  mm‘yS  ln  the  present  experiments.  p4  and  pi  were  adjusted 
accordingly  to  produce  shock  speeds  from  1100-1600  m/sec.  The  above 
lestnctions  were  not  prohibitive  since  they  provided  a postshock 
temperature  range  of  approximately  750-1400°K  and  densities  from 
about  0.35-G.75  the  density  of  air  at  STP  conditions. 

Another  postshock  parameter  which  had  to  be  considered  was  the 
vibrational  excitation  time  of  a gas  constituent  compared  to  the 
test  time  of  the  shock  tube.  The  excitation  time  had  to  be  shorter 
than  the  test  time  or  equilibrium  conditions  could  not  be  achie  -cd 
behind  the  wave.  The  test  time  ranged  from  100-300  usee,  therefore, 

air  could  be  investigated  for  lower  speed 
shocks.  At  the  higher  shock  velocities,  the  properties  of  N?  in 
air  were  investigated.  Pure  N2  could  not  be  examined  at  any  equilibrium 
conditions  in  the  present  tubo'because  of  the  long  excitation  iimes 
compared  with  the  test  time  restrictions. 
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IV.  RAMAN  INSTRUMENTATION  FOR  SHOCK  TUBE  DIAGNOSTICS 


A.  GENERAL  DISCUSSION 

The  application  of  Raman  spectroscopy  to  a shock  tube  environment 
is  in  general  a straightforward  process  of  obtaining  the  required 
scattering  data  at  the  proper  time  and  within  a short  time  interval. 
Also,  the  measurements  must  be  made  in  a small  region  within  the 
center  of  the  tube.  The  recording  time  must  be  short  in  order  that 
instantaneous  and  not  time  averaged  results  are  obtained  within  the 
test  time.  When  rate  measurements  are  being  made,  the  measurement 
time  must  also  be  short  compared  to  the  nonequilibrium  thermochemical 
times  associated  with  the  shock  heating  process.  The  scattering 
volume  must  be  small  to  avoid  boundary  layer  regions  of  the  flow 
which  would  complicate  the  calculation  of  postshock  gas  properties. 

A small  scattering  volume  dimension  parallel  to  the  shock  front  Is 
also  required  if  accurate  time  histories  are  to  be  obtained  behind 
the  wave  for  thermochemical  ratJ  measurements. 

The  particular  experimental  instrumentation  and  associated  optical 
apparatus  necessary  to  fulfill  the  above  requirements  are  discussed 
in  the  following  sections. 

B.  OVERALL  EXPERIMENTAL  SETUP 

A schematic  of  the  Raman  scattering  apparatus  for  the  present 
shock  tube  measurements  is  shown  in  Figure  5.  Figure  6 is  a photograph 
of  the  setup  which  emphasizes  the  optical  components  and  the  path 
of  the  laser  beam.  In  Figure  7,  a view  of  the  entire  experimental 
setup  is  shown  which  includes  the  spectrum  analyzer  and  recording 
oscilloscope. 

Since  the  shock  tube  generates  high  temperatures  and  pressures 
for  only  a few  hundred  microseconds  and  thermochemical  times  can  be 
orders  of  magnitude  less  than  that,  pulsed  lasers  are  best  suited 
for  Raman  measurements  in  these  transient  flow  facilities.  Hence, 
a 100  megawatt  pulsed  ruby  laser  with  a pulse  width  of  20  nanoseconds 
was  used  as  the  monochromatic  source. 

o 

The  ruby  radiation  at  6943  A traveled  along  a path  under  the  shock 
tube  and  in  a direction  parallel  to  the  tube's  axis.  The  beam  was 
then  rotated  90°  as  shown  in  Figure  5.  Two  right  angle  internal 
reflecting  prisms  were  actually  used  as  shown  in  Figure  6 to 
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properly  turn  the  horizontally  polarized  beam  and  also  rotate  the 
axis  of  polarization  so  that  it  was  in  the  plane  of  the  paper  as 
the  radiation  passed  through  the  test  volume.  The  1 centimeter 
diameter  beam  was  then  focused  with  a cylindrical  lens  into  the 
center  of  the  25  mm  diameter  tube.  A 6 mm  diameter  spatial  filter 
was  used  to  obtain  a 1 x 6 mm  beam  dimension.  Intersection  of  this 
beam  volume  with  the  0.5  x 6 mm  high  slit  image  of  the  spectrum 
analyzer  defined  a 0.5  x 6 x 6 mm  scattering  volume  shown  in  Figure  5. 

Thus,  the  measured  densities  and  temperatures  were  obtained  from  light 
scattered  within  the  above  volume.  This  small  spatial  resolution 
implied  that  good  time  resolution  was  obtained  behind  the  shock  wave 
and  also  insured  that  boundary  layer  effects  were  negligible. 

Comparison  could  then  be  made  of  Raman  data  with  simple  shock  tube  : , 

theory.  :•  ’ 

After  passing  through  the  test  region,  the  laser  radiation  was  T 

collected  in  a light  trap  to  avoid  detection  by  the  Raman  signal  : 

processing  instrumentation.  A small  portion  of  the  incident 
radiation  was  first  diverted  with  a beamsplitter  to  the  laser  power 
monitor  which  consisted  of  a photodiode  and  appropriate  attenuation 
filters.  The  power  monitor  was  required  in  order  to  record  pulse-  \ 

to-pulse  variation  in  the  laser  output.  Since  the  Raman  signal  is  J 

proportional  to  the  incident  laser  intensity,  this  variation  was 
ratioed  out  to  obtain  accurate  Raman  intensities.  Note  that  only  ] 

relative  intensity  measurements  are  required  so  absolute  laser 
power  measurements  are  not  necessary.  This  fact  will  be  discussed  f' 

in  more  detail  later  in  the  report.  \ . 

As  shown  in  Figure  8,  the  Raman  scattered  light  at  90°  to  the  | 

incident  laser  beam  was  collected  with  a f/2.5  lens  located  a focal  ; 

distance  of  13  mm  from  the  scattering  volume.  A f/5  lens  focused 
the  collected  radiation  on  to  the  slit  of  the  spectrum  analyzer. 

i 

C.  ELECTRO- OPTICA  L C CMP  CNENTS  I ! 

t j 
! 

The  scattered  light  collected  by  the  optics  must  be  analyzed  ; «; 

to  determine  its  spectral  content.  In  the  present  experiment  a ! 1 

3/4  meter  f/6.7  Spex  spectrometer  was  used.  A 600  gr/mm  grating 
blazed  at  lp  was  employed  in  conjunction  with  1 mm  slit.  This 
combination  resulted  in  an  apparatusofunction  which  was  nearly  tri- 
angular with  a 1/2  bandwidth  of  22. 5A  in  the  first  order.  The  true 
apparatus  function,  A(AA),  was  measured  with  a helium  neon  laser  and 
the  data  fit  with  mathematical  functions.  These  equations  were  coded 
into  the  Raman  theoretical  intensity  program  to  predict  the  unresolved 
spectra  as  a function  of  temperature  using  Eq.  (13). 

After  the  scattered  light  was  spectrally  analyzed  it  was  converted 
to  an  electrical  signal  by  an  RCA  31034A  photomultiplier  tube,  PMT, 
which  was  selected  for  a high  quantum  efficiency  of  19  percent  at 
8600A.  Before  reaching  the  PMT,  the  light  passed  through  a laser 
blocking  filter.  The  filter  used  in  the  present  study  was  a 3 mm 
thick  piece  of  RG-N9  black  glass  made  by  Schott  Optical  Glass,  Inc. 
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This  filter  transmitted  85  percent  of  the  incident  radiation  at 
wavelengths  greater  than  7700A  fi.e.,  the  Raman  region  of  interest) 
yet  blocked  ruby  laser  radiation  by  a measured  factor  of  nearly  10s. 

The  filter  was  used  to  eliminate  laser  radiation  which  was  not  removed 
by  the  spectrometer.  Removing  this  unwanted  light  was  necessary  to 
prevent  masking  of  the  weak  Raman  radiation  from  the  test  gas. 

Raman  signal  from  the  PMT  was  recorded  on  channel  #2  of  a Tektronix 
454  oscilloscope  using  an  EGfiC  Model  IT  100  inverting  transformer  to 
invert  the  negative  PMT  signal.  A SO  ft  feed-through  terminator 
followed  the  inverting  transformer  to  provide  fast  response  times  in 
order  to  accurately  record  the  20  nanosecond  pulses.  A positive  output 
signal  from  the  Korad  KD-1  photodiode  incident  laser  intensity  monitor 
was  also  terminated  in  50  ft  and  recorded  on  the  same  oscilloscope  by 
feeding  this  signal  into  channel  #1  and  putting  the  scope  in  the  "add" 
mode.  Since  there  was  approximately  a 100  nanosecond  transit  time  of 
the  electrons  through  the  dynode  channel  of  the  PMT,  then  this  signal 
arrived  after  the  photodiode  signal.  Hence,  both  were  recorded  on  the 
same  instrument.  An  example  oscilloscope  trace  recorded  in  this  fashion 
is  shown  in  Figure  9.  A trace  of  this  type  was  obtained  for  each  laser 
firing.  Data  points  then  correspond  to  ratioing  the  Raman  signal, 

Ir,  to  the  laser  signal,  II,  taken  off  the  same  trace. 

The  oscilloscope  sweep  was  triggered  by  utilizing  the  negative 
energy  pulse  which  is  also  provided  by  the  KD-1  photodiode  circuit  at 
the  time  the  laser  is  fired. 

D,  COINCIDENCE  TIMING 

As  mentioned  above,  the  laser  had  to  be  fired  at  the  proper  time 
behind  the  shock  front.  At  first  the  solution  to  this  problem  seemed 
to  be  a simple  matter  of  using  a single  transducer  to  provide  an 
output  when  the  shock  wave  passed  a certain  point.  This  output  could 
be  suitably  delayed  to  fire  the  laser  at  the  proper  time.  However, 
the  nature  of  a ruby  laser  firing  sequence  complicated  matters  somewhat. 
The  ruby  rod  must  be  optically  pumped  with  a flash  lamp  for  approximately 
1 millisecond  before  a giant  pulse  can  be  obtained  by  Q-switching  the 
cavity.  Hence,  using  a single  trigger  meant  that  flash  lamp  triggering 
would  occur  when  the  shock  wave  was  only  about  half  way  down  the  tube. 
Small  variations  in  shock  speed  were  amplified  over  a long  distance 
and  a 50  psec  jitter  was  encountered  in  the  time  that  the  laser  was 
fired  with  respect  to  passage  of  the  shock  front.  The  actual  position 
of  the  wave  with  respect  to  the  laser  firing  was  easily  obtained  to  a 
high  degree  of  accuracy  by  using  the  KD-1  photodiodes’  negative  energy 
pulse  to  start  a timer  which  was  stopped  by  passage  of  the  wave  by 
position  #3  in  Figure  5.  However,  using  a single  transducer,  the 
actual  measured  time  would  be  a large  random  scattered  value  around  a 
preset  time  rather  than  the  value  picked  by  the  experimenter. 

To  remedy  this  situation,  one  transducer  located  near  the  diaphragm 
(i.e.,  transducer  DT  in  Figure  10)  and  delay  was  used  to  fire  the  lamp 
while  a second  transducer  and  delay  fired  the  Q-switch.  This  second 


I 


26 


wmmmmm 


NSWC/WOL/TR  75-53 


transducer  (i.e.,  # 1 in  Figure  10)  was  located  near  the  scattering 
volume  in  order  to  obtain  a jitter  r.f  only  ±2  microseconds.  Other 
details  of  the  timing  and  coincidence  circuitry  are  shown  in 
Figure  10.  Each  trigger  circuit  consisted  of  a pressure  transducer 
and  a cathode  follower  amplifier  which  was  used  to  drive  a trigger 
generator.  The  5 volt  output  of  the  trigger  generator  started  a preset 
delay.  The  15  volt  output  of  the  delay  provided  a trigger  signal 
compatible  with  the  laser  circuitry.  Using  this  two  trigger  method, 
the  transducer  signals  were  suitably  delayed  to  obtain  both  maximum 
pumping  of  the  laser  rod  for  efficient  lasing  and  proper  orientation 
of  the  shock  front  with  respect  to  the  viewing  port.  The  Q-switch 
transducer  was  also  used  to  start  two  Monsanto  counter  timers  for 
shock  speed  measurements  as  mentioned  above. 
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V.  SHOCK  TUBE  RAMAN  DATA 


A.  TIME  HISTORY  OF  SCATTERING  INTENSITY 

Raman  scattering  intensities  were  monitored  as  a function  of  time 
as  a shock  wave  passed  the  viewing  port-  Incident  laser  intensities 
were  also  recorded  and  Ip/I],  ratios  were  obtained.  Results  of  these 
measurements  are  shown  in  Figure  11-  Each  data  point  corresponded  to 
a single  firing  of  the  shock  tube.  The  laser  was  also  fired  once  with 
each  shot  but  with  a variable  delay.  0Hence,  the  first  portion  of  the 
data  (i.e.,  open  circles)  at  A = 8280A  in  Figure  11  was  scattering 
from  nitrogen  ahead  of  the  wave  at  the  initial  driven  gas  conditions 
(i.e.,  P = P^)*  The  final  data  points  correspond  to  scattering 
from  shock  heated  nitrogen.  At  62  usee  before  the  wave  passes  position 
$3  (see  Figure  5)  , the  Raman  scattering  intensity  increases  abruptly. 
This  increase  represents  the  passage  of  the  wave  through  the  scattering 
volume.  The  constant  intensity  conditions  behind  the  wave  lasted  for 
approximately  150  usee.  After  that  time  a second  abrupt  change  in 
scattering  intensity  occurred.  This  second  increase  was  most  likely 
due  to  passage  of  the  contact  surface  through  the  scattering  volume. 

The  measured  test  time  compared  quite  closely  to  the  estimated  value 
for  the  experimental  conditions  listed  in  the  figure. 

To  insure  that  the  output  in  Figure  11  was  not  some  false  signal, 
the  intensity  at  A = 8310A  was  also  measured  as  a function  of  time 
(i.e.,  square  symbols  in  Figure  11).  Note  that  negligible  output  was 
obtained  at  this  wavelength  which  was  just  outside  the  N!2  Raman  region. 
Hence,  the  recorded  output  at  8280A  was  truly  due  to  Raman  scattering 
from  nitrogen. 

Also  shown  in  Figure  11  (i.e.,  diamond  symbols)  is  the  amount  of 
scattering  obtained  from  600  torr  of  nitrogen  at  static  conditions. 

This  pre-shot  data  was  used  as  a calibration  bench  mark  to  determine 
the  absolute  density  prior  to  and  behind  the  wave.  The  data  in 
Figure  11  was  reduced  to  true  density  ratios  by  dividing  each  data 
point  by  the  average  intensity  value  of  Ir/I^  = *13  ahead  of  wave. 

The  data  was  also  corrected  for  temperature  effects  as  discussed 
in  Section  1 1 - D . This  temperature  correction  will  be  further  clarified 
in  the  following  sections.  The  measurements  in  Figure  11  were  also 
adjusted  by  subtracting  62  nsec  from  the  recorded  times.  The  results 
of  this  entire  data  reduction  process  are  shown  in  Figure  12.  Also 
shown  in  the  figure  is  the  theoretical  plot  for  density  ratio  based 
on  constant  ratio  of  specific  heat  of  y = 1.4.  This  theory  was  picked 
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because,  for  pure  nitrogen  under  the  present  experimental  conditions, 
the  vibrational  mode  would  not  be  excited  in  the  short  test  time 
encountered  in  this  experiment.  Note  that  the  measured  ratio  was 
10  percent  greater  than  the  theoretical  prediction  for  this  particular 
case.  This  disagreement  was  reasonable  considering  that  the  data  in  ’ 

Figure  12  were,  to  the  author's  knowledge,  the  first  Raman  measurements 
behind  a shock  wave  and  much  was  learned  about  improving  the  experimental 
procedure  following  this  initial  data  gathering  sequence. 

B.  DETERMINATION  OF  VIBRATIONAL  EXCITATION  TIME  FOR  OXYGEN  IN  AIR 

I 

As  discussed  in  Section  II-E,  the  vibrational  excitation  time  for 
a gas  can  be  determined  by  examining  the  Raman  scattering  profile  as 
a function  of  time  behind  the  wave.  This  time  history  data  was 
obtained  for  oxygen  vibrational  excitation  in  air  at  a shock  velocity 
of  approximately  1100  meters/sec  and  initial  air  pressure  of  150  torr. 

The  data  consisted  of  measuring  the  ratio  of  the  Raman  scattering 

intensity  at  Aj  = 7762.A  to  the  incident  laser  intensity  at  various 

times  relative  to  the  passage  of  the  shock  front.  The  results  are 

shown  in  Figure  13.  Also  plotted  for  illustrative  purposes  in  the 

figure  are  the  constant  density  profile  shapes  for  various  shock  1 

conditions  and  the  positions  on  the  time  history  plot  where  these  ; 

conditions  prevail.  Note  that  the  first  profile  represents  static 

gas  at  ambient  conditions  ahead  of  the  wave.  The  Raman  profile  for 

this  condition  is  symmetric  indicating  negligible  upper  level  rotational 

or  vibrational  excitation.  The  second  condition  labeled  (21  on  the 

plots  is  indicative  of  the  state  of  the  gas  immediately  behind  the 

shock  wave.  Here  the  rotational  temperature  has  increased  to  at  least 

the  equilibrium  value  behind  the  wave  while  the  vibrational  temperature 

still  remains  at  the  low  ambient  value.  The  increase  in  rotational 

temperature  results  in  a change  in  the  profile  to  a more  asymmetric 

shape  as  shown  in  the  figure.  The  asymmetry  is  caused  by  the  increased 

population  of  the  upper  rotational  states  of  the  oxygen  molecule. 

Note  that  the  increase  in  Raman  intensity  between  (1)  and  (2)  appears 
to  be  only  a factor  of  3 due  to  the  profile  changes  alone.  However, 
the  data  indicates  a change  closer  to  a factor  of  12.  The  apparent 
difference  between  these  two  factors  is  explained  by  the  fact  that  the 
density  also  increases  behind  the  wave  by  a factor  of  4, 

Finally,  the  vibrational  mode  excitation  increases  with  time  behind 
the  wave  until  the  third  condition  is  reaches.  This  final  state 
represents  fully  excited  vibration  and  rotation  which  results  in  a 
further  asymmetry  m the  Raman  profile. 

This  final  Raman  intensity  increase  at  A]  was  approximately  a 
factor  of  2 as  noted  from  the  theory  and  experimental  data.  The  time 
between  condition  (2)  and  (3)  is  defined  as  the  vibrational  excitation 
time  for  oxygen  in  air. 

Also  shown  in  Figure  13  is  the  theoretical  Raman  intensity  ratio 
represented  by  the  solid  line.  Note  that  a dashed  line  was  used  between 
(2)  and  (3)  to  indicate  that  this  path  is  not  being  exactly  represented. 
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The  prediction  would  be  quite  complex  for  several  reasons.  There 
would  actually  be  a rotational  temperature  overshoot  occurring  just 
behind  the  wave  and  the  exact  vibrational  temperature  and  rotational 
temperature  would  be  difficult  to  predict.  In  actuality,  there  is  no 
real  heed  to  Know* the  exact  path  because  only  the  end  points  are 
required  for  this  excitation  time  measurement.  In  Figure  13  the 
theory  has  been  implicitly  corrected  for  temperature  effects  by 
accounting  for  profile  changes  with  temperature,  in  fact,  the 
existence  of  temperature  effects  is  the  reason  behind  the  ability  to 
measure  the  vibrational  excitation  process. 

Note  that  the  measured  Raman  intensities  at  the  end  points  agree 
quite  well  with  the  predicted  values.  However,  the  measured  time 
of  36  psec  shown  in  figure  13  is  30  percent  greater  than  the  theoretically 
predicted  value  of  27  nsec. 


The  theoretical  excitation  time  was  Computed  utilizing  experimental 
results  determined  in  other  shock-tube  experiments  using  interferometric 
techniques . First  the  times  had  to  be  converted  to  lab  fixed  coordinates 
from  wave  fixed  values  since  the  Raman  measurement  was  in  the  lab 
fixed  coordinate  system.  To  convert  from  one  system  to  another  the 
f'dlipwihg  fPrmuia  is  used: 


pT 


w 


Lab  p^Cp^/pp 


(17) 


where  t 
shock  p 


and  xr  are  the  lab  and  Wave  fixed  values,  p?  is  the  post- 
sure  and  (p2/Pl)  is  the  density  ratio  across“the  shock. 


• w for  oxygen  excitation  in  air  (i  .e. , tq  _^jr)  * 5 8*vcn  by: 


pT02-Air 


(18) 


In  Eq.  (18).  pj  is  the  mole  fraction  of  gas  i in  air.  pr0  and 
prn  v are  given  by:  u2"  2 

u 2 'n2 


PTn  n . 1056  CT'1/S-.03)-8 
-2  ' u2 


(atm- sec] 


(19) 


P1n.  -n.  ~ 1 • 4 pip  h 


(20) 


[atm- sec] 
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Eq,  (19)  was  obtained  from  shook  tube  data  .by  White  and  Millikan.15  j| 

Eq,  (20)  was  obtained  from  observations  of  Blackman16  showing  that  jf 

N*2  is  40  percent  as  effective  as  ©2  in  the  V-T  excitation  of  0? . § 

C.  TEMPERATURE  AND  DENSITY  DETERMINATION  USING  SPECTRAL  SCAN  I 1 

METHOD 

As  discussed  in  Sections  1 1 - D and  II-E,  the  thermochemical  state 
of  a gas  can  be  determined  by  measuring  its  Raman  scattering  spectrum. 

This  intensity  versus  wavelength  data  can  then  be  used  to  obtain  the 
temperature  and  density  of  various  species  in  a gas  mixture.  _ In  this 
section,  the  postshock  temperature  and  density  measurements  of  Cb 
and  NS  in  air  are  discussed. 


The  objective  was  to  determine  equilibrium  gas  properties  and 
compare  the  Raman  measurements  to  equilibrium  theoretical  calculations. 

To  insure  that  indeed  local  thermodynamic  equilibrium  (LTE)  did  exist 
during  the  test  time,  calculations  were  made  of'the  vibrational 
excitation  times  and  these  computed  values  were  compared  to  the  expected 
test  time  in  this  particular  experiment.  As  discussed  in  Section 
1 1 1 - D , and  V-A , pure  NS  could  not  be  examined  under  eq,-lli  hrium  conditions 
due  to  its  slow  excitation  time.  However,  calculations  showed  that  O2 
in  air  was  sufficiently  rapid  such  that  equilibrium  conditions  could 
be  obtained  over  a wide  range  of  conditions.  Also,  H2  in  a; r could 
be  examined  at  higher  shock  velocities  due  to  the  enhanced  late  of 
the  N2  V-T  process  in  the  presence  of  oxygen.  For  LT^  conditions, 
the  temperatures  of  the  various  modes  are  equal  to  the  translational 
temperature  of  the  gas  (i.e.,  Tr  = Tv  = T ) . As  discussed  at  the  end 
of  Section  II-C,  a broad  apparatus  function  can  be  employed  in  this 
case*  The  following  data  was  thus  obtained  uging  the  measured  apparatus 
function  xvith  a width  at  half  height  of  22. 5A. 


To  obtain  accurate  Raman  temperature  and  density  measurements  using 
LRS  in  Shock  tube  experiments,  various  precautions  had  to  be  observed 
during  the  data  taking  sequence.  The  first  step  was  to  make  sure  that 
the  scattering  volume  was  properly  aligned  with  the  spectrometer  slit, 
this  was  accomplished  by  translating  the  f/2.5  lens  in  Figure  S in  a 
direction  perpendicular  to  the  slit  height.  This  adjustment  traversed 
the  narrow  1 mm  beam  image  across  the  1 mm  slit  width.  Raman  intensity 
at  a fixed  wavelength  and  pressure  was  then  obtained  as  a function  of 
lens  position  prior  to  a data  gathering  sequence.  The  position  of  peak 
scattering  intensity  was  found  and  the  lens  location  set  at  that  position. 
The  intensity  versus  position  profile  was  relatively  flat  near  the 
peak  value.  However,  periodic  cheeks  showed  that  small  adjustments 
were  required  to  maintain  proper  alignment.  The  stiffener  which  was 
added  to  the  tube  helped  to  reduce  misalignment  problems. 


^White,  D,  R.  and  Millikan,  R.  C.,  "Vibrational  Relaxation  in 
Air,"  AIAA  Journal  2,  pp.  1844-1846,  October  1964 

^Blackman,  V,,  "Vibrational  Relaxation  in  Oxygen  and  Nitrogen," 
■J.  Fluid  Mech.  1,  pp.  61-85,  May  1956 
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The  next  step  was  to  obtain  a Raman  spectral  scan  at  static 
conditions  for  the  particular  gas  to  be  measured  in  the  postshock 
state.  This  data  consisted  of  ratioing  Raman  scattering  intensities 
to  incident  laser  intensity  values  (i.e.,  Tg/I]J  for  various 
wavelengths  throughout  the  Raman  region.  The  results  of  such  a scan 
for  nitrogen  in  air  are  shown  on  the  upper  plot  in  Figure  14.  This 
raw  data  was  obtained  by  evacuating  the  tube  to  a pressure  of  200  torr 
of  air  prior  to  a given  testing  period.  Each  data  point  (i.e.,  open 
symbol)  corresponds  to  a single  laser  firing  but  at  a different 
indicated  spectrometer  wavelength  position,  Aj.  Note  that  some  residual 
intensity  not  associated  with  Raman  scattering  was  present  outside  the 
N!2  Raman*  region  (i.e.,  iBp  = .2).  This  background  light  was  probably, 
due  to  photons  produced  from  fluorescence  of  the  liquid  absorber  used 
in  the  laser  dump.  The  measured  Raman  intensities  had  to  be  corrected 
for  this  effect  by  subtracting  Igp  from  the  data  in  Figure  14  to 
obtain  1^,  the  true  calibration  intensity.  After  this  background  light 
subtraction  was  accomplished,  the  data  was  normalized  -with  the  peak 
intensity  value,.  ImaX£»  using  the  following  formula: 


i 


I 


(VIP-1*, 


I 1 

maxj,  max^ 


1 

(21)  1 


Note,  the  value  of  tmaxc  Eq*  (21)  was  obtained  by  drawing  a smooth 
line  through  the  .data  points  (i.e.,  in  -Figure  ,14 , Imax^  = 2 .30- . 20=2 . 10) . 
The  purpose  of  obtaining  the  normalized  static  scan  may  not  be 
obvious  at  this  time.  However,  to  this  point  in  the  data  gathering 
sequence,  the  discussion  has  centered  around  indicated  spectrometer 
wavelength  positions,  Aj.  These  values  are  not  a priori  equal  to  the 
actual  wavelengths  appearing  at  the  exit  slit  of  the  spectrometer. 

In  fact,  it  would  be  fortuitous  if  they  were,  A wavelength  calibration 
procedure  using  standard  lamp  sources  was  helpful  in  determining  an 
approximate  value  of  the  wavelength  correction,  A,  given  by: 


A =Ai‘-Aa  (22) 

Here  A ^ is  the  actual  wavelength  position  while  Aj  is  the  position  that 
is  indicated  on  the  wavelength  dial  of  the  instrument.  But,  A is  a 
small  but  non-negligible  function  of  A^  itself  so  the  calibration  source 
wavelength  would  have  to  be  located  almost  identical  to  the  Raman 
wavelength.  This  was  not  practical  for  the  present  experiment,  so, 
instead,  the  Raman  data  from  Eq.  (21)  at  known  static  conditons  was 
used  to  accurately  determine  the  Aoto  be  used  to  correct  the  spectro- 
meter wavelength  settings.  A = 20A  was  determined  from  the  normalized 
form  of  the  data  in  Figure  14.  For  purposes  of  comparing  pre-  and 
postshock  data,  the  static  Raman  output  in  Figure  14  was  further 
reduced  to  values  indicative  of  initial  preshock  conditions 


*“*  O 
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(i.e.,  p - pp.  This  was  accomplished  by  multiplying  the  intensities 
that  were  corrected  for  background  light  by  a ratio  of  the  initial 
pressure,  pi,  to  the  pressure  at  which  the  static  scan  data  was  obtained 
(i.e.,  p = Pcalib) * Hence: 


1 1M  = TC  X 


x - -PL. 
C Pcalib. 


Where  I are  the  measured  static  Raman  scan  values  referenced  to  the 
initial  preshock  conditions.  The  results  of  this  data  reduction  are 
shown  as  open  circles  in  Figure  15. 

To  compare  this  static  scan  data  to  a theory,  Raman  spectral0 
intensities  were  computed  using  Eq.  (13)  at  T = 296°K.  The  22. 5A 
experimentally  determined  apparatus  function  was  used  in  these 
calculations.  The  computer  output  intensities,  T are  normalized 
to  give  a peaik  intensity  of  1.,  These  values  were  then  adjusted  to 
the  preshock  experimental  conditions  using: 


t1T  = t-v  (T  = 296°K)  x 


calib 


The  theoretical  Raman  profile  at  the  initial  static  preshock 
conditions,  is  the  solid  curve  plotted  in  Figure  15.  Ii-p  was 

plotted  in  this  figure  to  illustrate  that  the  computer  calculation 
compared  well  with  room  temperature  data  and  also  to  show  that  the 
value  of  A is  reasonably  accurate. 

Having  performed  the  necessary  alignment  and  wavelength  calibration; 
the  raw  shock  tube  data  was  then  obtained.  The  procedure  was  the  same 
as  for  the  static- scan  data  except  that  the  passage  of  the  shock  front 
triggered  the  laser  to  fire  at  the  proper  time  behind  the  wave.  In 
the  present  experiment  this  time  was  taken  as  the  value  after  the 
measured  excitation  time  and  before  the  useful  test  time.  The 
results  are  shown  on  the  lower  plot  in  Figure  14  for  shock  heated 
nitrogen  in  air.  The  postshock  scan  data  must  again  be  corrected  for 
background  light  by  subtracting  Ig*.  = .27  from  the  measured  Raman  to 
incident  laser  intensity  ratios  in'rigure  14  to  obtain  Ig,  the  tube 
scan  intensity.  Tg  versus  wavelength  data  is  then  plotted  in  Figure  15 
using  open  square  symbols.  The  postshock  measurements  can  now  be 
compared  to  static  data  at  preshock  conditions. 

To  relate  this  scan  data  to  actual  density  and  thus  to  obtain 
density  ratios,  an  accurate  standard  calibration  is  required  (i^e. , 
more  accurate  than  the  static  scan  data  previously  obtained).  Such 
a calibration  simply  involves  measuring  the  Raman  scattering  intensity 
at  a knwon  pressure  (i.e.,  p * Pgjp)  at  ambient  temperature  conditions. 
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Furthermore  this  standard  intensity,  Ig^,  should  be  measured  at  a 
fixed  wavelength  location  (i.e.,  = AgTp) . 

At  this  point  it  would  be  appropriate  to  say  something  about  the 
magnitude  and  meaning  of  the  Ir/I^  values  found  in  Figure  14.  Note 
that  in  Figure  9,  Ir/I^  ~ 2.27  was  obtained.  The  point  to  be  made 
here  is  that  this  number  is  completely  arbitrary  and  the  absolute 
ratio  need  not  be  determined  as  long  as  nothing  unknowingly  happens 
to  alter  the  sensitivity  of  the  intensity  monitoring  equipment  between 
the  shock  tube  tests  and  the  standard  calibration  sequence.  To  make 
sure  the  sensitivity  is  not  changing  with  time,  the  standard 
calibration  test  should  be  made  periodically  between  sets  of  shock  tube 
firings.  If  this  relative  Raman  calibration  output  (i.e.,  Igyp)  at  a 
known  ambient  condition  is  constant  with  time  during  the  shock  heated 
tests,  then  the  sensitivities  have  remained  constant  and  the  calibration 
intensities  can  then  be  used  to  obtain  exact  results  from  the  arbitrary 
intensity  values  shown  in  Figures  14  and  15 ; The  average  calibration 
intensity  value  obtained  prior  to  and  during  the  shock  tube  data 
sequence  is  also  plotted  as  2.18  in  Figure  14  (i.e.,  solid  circular 
symbol)  along  with  the  standard  deviation  bars  representing  scatter 
about  this  value.  Again  this  calibration  must  be  corrected  for 
background  light  intensity  by  subtracting  Irc  from  the  average  value 
to  obtain  the  actual  standard  calibration  intensity  (i^e.,  TsTD  ~ 

2 . 18- . 20=1 . 98) . Also  note  from  Figure  14  that  the  average  intensity 
value  fits  quite  well  with  the  raw  scan  data  obtained  previously. 

I STD  was  then  altered  to  reflect  preshock  conditions  so  it  could  be 
used  to  determine  the  density  ratio  across  the  shock.  Hence,  the 
standard  intensity  at  state  (T)  conditions  was  given  by: 


I 


1 


STD 


TSTD 


x 


Pi 

PSTD 


(25) 


(i.e.,  I, 

»STD 


where:  pg,^ 


1.98  x 1/4  = .495  from  Figure  14) 


Pcalib 


200  torr 


Note  that  Eq.  (25)  represents  the  same  correction  procedure  that  was 
performed  on  the  scan  data  with  Eq.  (23).  llg-rn  = -495  is  then 
plotted  in  Figure  15  for  comparison  to  the  static  scan  data  and  also 
to  reduced  values  of  postshock  data.  The  standard  wavelength,  A 
was  828  2 A . 6 iV 

The  postshock  Raman  output  in  Figure  14  was  then  curve  fit  in  order 
to  obtain  the  temperature  and  density  -of  nitrogen  in  the  shock  heated 
air.  The  trial  curve  fits  produced  residual  error  sums  whose  magnitude 
is  based  on  how  well  the  theory  agreed  with  measured  data.  Two 
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variables  were  adjusted  to  obtain  a minimum  residual  value.  These 
variables  were  gas  temperature  and  maximum  intensity,  I.MAXg*  of  che 
theoretical  profile.  l,MAXs  Kas  used  to  normalize  the  scan1  data  to 
1 in  order  to  compare  the  measured  values  to  the  Raman  theory  which 
was  also  normalized  to  1.  Hence,  Raman  theory  refers  to  results  of 
Eq.  (13)  normalized  to  the  peak  intensity  value.  The  curve  fit 
results  give  postshock  temperature,  hence  the  temperature  ratio, 
directly.  However,  to  obtain  the  shock  density  ratio,  further  data 
reduction  is  required  using  the  iMAXg  value  obtained  from  the  curve 
fit  process.  This  further  data  reduction  is  straightforward  and 
consists  of  finding  the  best  fit  value  for  the  postshock  intensity 
at  ^STD  = 8282A.  To  accomplish  this  the  normalized  theoretical  Raman 
spectral  profile  which  best  fit  the  Raman  data  was  multiplied  by 
IiMAXcj  to  obtain  Iy,  where: 


tbest  fit^ 


Results  of  the  calculations  using  Eq.  (26)  are  then  plotted  for 
comparison0 with  the  postshock  data  in  Figure  15.  Also,  the  value  of 
!-r  at  8282A  was  found  to  be  1.89  and  is  labeled  in  the  figure. 

At  long  last,  the  density  ratio  can  be  computed,  without  correcting 
for  temperature  effects,  by  ratioihg  the  two  intensity  values  at 
AsTI)  in  Figure  15.  Hence: 


uncorrected 


If (A  Astd) 
* 1 

\STD 


(i.e.,  P7/p^  = 1.89/. 495  = 3.82  from  data  in  the  figure) 

Finally,  the  true  density  ratio  can  be  determined  from  the  reduced 
data  by  correcting  the  above  for  temperature  effects.  The  reason  that 
the  temperature  was  hot  accounted  for  in  the  previous  procedures  is 
that  the  theoretical  computations  used  above  were  always  normalized 
to  a peak  intensity  Of  1 for  convenience  sake  in  performing  the  curve 
fit  process.  However,  as  Figure  3 depicts,  the  peak  intensity  (and 
hence  all  other  intensity  values)  varies  as  a function  of  temperature 
at  constant  density.  Hence,  this  temperature  effect  must  be  accounted 
for  in  the  data  reduction  process.  This  was  accomplished  by  multiplying 
the  uncorrected0value  by  a ratio  of  the  theoretical  Raman  output  at 
A = Ag'j’p  = 8282A  for  room  temperature  fi.e.,  calibration  temperature) 
conditions  to  the  theoretical  output  at  the  postshock  temperature. 

Hence : 


uncorrected 


Ir(T=T 


BEST  FIT' 
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The  theoretical  Raman  intensity  as  a function  of  temperature, 

Ir(A).=  in  Eq.  (13j,  is  plotted  in  Figure  16  for  A = Ap  = 8282A. 

From  the  figure,  at  T = Tgggy  pjy.  Ip  = >743  and  at  T = 29b°K,  Tp  = 1. 
Using  Eq.  (28),  the  temperature  correction  follows  directly  (i.e., 

\ / o?  = 3.82  x 1/.743  = 5.14)*  This  value  is  26  percent  higher  than 
the  uncorrected  density  ratio  which  indicates  the  importance  of  the 
temperature  correction  process. 

The  above  data  reduction  procedure  seems  laborious  and  complex. 

Its  description  could  have  been  simplified  somewhat  but  a lot  of  the 
details  and  some  of  the  physical  meaning  of  the  data  would  have  been 
lost.  Basically,  Figure  14  represents  the  raw  data  containing  arbitrary 
but  consistent  output  sensitivities  between  calibration  and  test 
conditions.  The  raw  data  was  corrected  for  background  light  and  then 
referred  to  the  initial  conditions  prior  to  the  shock  wave  (i.e., 
p = Pj) . In  Figure  15,  the  pro-  and  postshock  arbitrary  intensity 
data  were  compared  on  a single  plot.  Temperature  was  obtained  by 
a simple  curve  fit  of  the  arbitrary  intensity  values.  The  density 
determination  was  more  involved.  Accurate  pre-  and  postshock  intensities 
at  A = Agyp  were  ratioed  to  obtain  the  uncorrected  density  ratio. 

This  ratioing  of  arbitrary  intensities  results  in  a real  density 
ratio  and  hence  actual  postshock  density  simply  because  the  preshock 
arbitrary  intensity. corresponded  to  a known  density  value.  The  necessary 
factors  for  correcting  the  density  ratio  for  temperature  effects  were 
obtained  from  Figure  16. 

D.  TEMPERATURE  AND  DENSITY  MEASUREMENT  USING  2 -COLOR  TECHNIQUE 

In  many  ways,  the  2-color  (i.e.,  2 -wavelength)  technique  is  similar 
to  the  spectral  scan  method  in  that  the  2-color  data  essentially 
consists  of  only, a three  point  wavelength  scan.  Two  intensity  values 
are  obtained  within  the  Raman  spectral  region  and  a third  point  at 
a wavelength  outside  the  Raman  band  in  order  to  determine  background 
light  intensity.  The  idea  was  to  obtain  more  accurate  data  at  the 
fewer  number  of  points  by  measuring  Raman  and  background  intensities 
at  each  wavelength  for  three  to  five  shock  tube- laser  firings.  The 
rationale  behind  a 2-color  technique  is  that  the  spectrometer  could 
be  replaced  with  a series  of  filters  and  photo-detectors.  This  would 
improve  the  signal  gathering  capabilities  of  the  optics  since  filter 
systems  can  have  a very  low  f /no.  And,  using  three  filters,  all  the 
necessary  experimental  data  to  determine  density  and  temperature 
could  be  obtained  with  a single  laser  firing. 

The  same  precautions  regarding  alignment  of  the  scattering  volume 
with  the  spectrometer  slit  were  required  for  the  2-color  data  gathering 
sequence.  Also,  a Raman  spectral  scan  at  static  conditions  was 
again  necessary  to  convert  spectrometer  wavelengths  to  true  wavelength 
-values.  The  results  of  this  scan  arc  shown  in  the  upper  portion  of 
Figure  17.  Note  that,  the  peak  intensity  value  of  2.51  is  somewhat 
different  than  in  Figure  14  indicating  a small  change  in  the 
sensitivity  of  the  system  in  the  time  between  the  two  types  of  data 
gathering  sequences-.  T-be  results  of.  Figure  17  were  corrected  for 


background  light  intensity  and  normalized  using  F,q.  (21),  where 
UiAXp  = ‘ 2 , 51  - . 20=2 . 31 . The  value  of  A determined  from  this  static 
scair'procedure  was  2oX,  which  was  identical  to  the  previous 
wavelength  correction  factor. 

As  before,  I was  obtained  using  Eq . (23)  and  the  results  of  this 
data  reduction  procedure  are  shown  as  open  circles  in  Figure  18. 

Also  plotted  as  the  solid  line  through  this  data  is  the  computed 
intensity  profile  at  the  calibration  temperature  and  pressure 
conditions  using  Eq.  (24). 

The  standard  calibration  results  taken  prior  to  and  during  the 
2-color  data  sequence  are  also  shown  in  Figure  17  with  a solid 
circular  symbol.  Using  Eq . (25),  the  standard  intensity  was  referred 
to  postshock  conditions  (i.e.,  Iic-pp  = (2. 28-. 20)  x 1/4  = .520). 
llcjTD  also  plotted  in  Figure  18  again  using  a solid  circular  symbol. 

The  shock  tube  2-color  data  was  then  obtained  and  the  three 
intensity  values  are  shown  plotted  as  Ip  = 1.98,  1 i = 1.33  and 
Ig  = .43  with  square  symbols  in  the  lower  portion  of  Figure  17.  The 
standard  deviation  of  the  data  at  these  wavelength  values  is  also  shown  i 1 
in  the  figure  using  error  bars.  The  smooth  curve  drawn  through  the  J 

postshock  data  points  in  Figure  17  is  the  expected  Raman  output  and 
is  presented  for  illustrative  purposes  only  since  no  curve  fitting  is 
required  for  the  2-color  technique.  As  the  profile  depicts,  X]  was  j 

picked  near  the  peak  of  the  Raman  profile  while  A2  was  located  in  the 
blue  shifted  region  of  the  Raman  spectrum.  The  choice  of  these  : , 

wavelength  positions  is  not  completely  arbitrary.  They  were  picked  so  j 
that  the  intensity  at  Ap  would  be  somewhat  insensitive  to  temperature 
effects  (i.e.,  Ijp(Ap)  ~ const,  x density).  Further,  A2  was  chosen  so 
that  the  ratio  Ir(A2)/Ir(Ai)  would  be  most  sensitive  to  temperature  ; 

changes. 

From  this  point,  the  data  reduction  method  departs  rather  markedly 
from  that  of  the  spectral  scan  technique.  The  measured  postshock 
Raman  intensities  are  again  corrected  for  background  light  by 
subtracting  Ig  from  Ip  and  I2  to  obtain: 


W ' h ' h - ^ 

1 r(A 7)  - I7  - Ip  = .90 


(29) 


The  postshock  2-color  intensity  values  in  Eq . (29)  arc  plotted  in 
Figure  18  as  solid  square  symbols.  These  two  values  (i.e.,  2-color 
intensities)  in  conjunction  with  Ilgyp  = .523  arc  all  the  reduced  data 
needed  to  compute  temperature  and  density  using  the  2-color  method. 
Note,  that  for  convenience,  Ap  = A^-pp.  The  solid  curve  drawn  through 
the  postshock  data  in  Figure  18  represents  the  theoretical  Raman 


I 

pi 
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profile  for  the  measured  postshock  temperature  conditions  and, 
again,  it  is  only  shown  for  illustrative  purposes. 

To  obtain  temperature  from  the  above  2-color  data,  the  values 
in  Eq.  (29)  are  ratioed  to  give: 


1 r( A 2^ 


RA  1 • 


1.55 


= .58 


In  Figure  19,  the  theoretical  value  of  this  ratio  is  computed  for0the 
particular  wavelengths  used  in  the  present  study  (i.e.,  Ai  = 8285A, 

A2  = 8260A) . Finding  the  value  of  the  measured  ratio  of  .58  from 
Eq.  (30)  on  the  vertical  axis  of  Figure  19  and  moving  horizontally 
until  the  18260/^8285  curve  is  intersected,  one  determines  directly 
that  the  measured  temperature  is  1440°K. 

Further,  the  uncorrected  density  ratio  is  given  by  an  equation 
similar  to  Eq . (27)  as: 


ft). 


uncorrected 


I CAX) 
iSTD 


(i.e.,  p2 /pj  = 1.55/. 52  = 3.98  from  data  in  Figure  18) 

Again,  the  true  density  can  be  determined  by  correcting  for 
temperature  effects.  This  was  accomplished  by  multiplying  the 
uncorrected  valug  by  a ratio  of  the  theoretical  Raman  output  at 
Ai  = A5'pj)  - 8285a  for  standard  ambient  temperature  conditions  to 
the  theoretical  output  at  the  measured  postshock  temperature.  Hence 


i -ft). 


uncorrected 


TR^T~TROaP 
1 R ^T=^MEASURErP 


This  process  is  analogous  to  that  of  Eq.  (28).  The0theoretical  Raman 
intensity  as  a function  of  temperature  at  A2  - 8285A  is  also  plotted 
in  Figure  19.  From  the  figure,  at  T= Treasured*  Ir  = -615  and  at 
T = 295°K,  IR  = .913.  Using  Eq . (32),  tne  temperature  correction 
follows  directly  (i.e.,  p7/p^  = 2.98  x .913/. 615  = 4.42) 
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VT.  DISCUSSION  OF  RESULTS 


Excitation  of  pure  O2  and  .\’2,  as  well  as  S2  and  O2  in  air  have  J 

been  discussed  in  this  report.  To  clarify  statements  made  concerning 
these  rates,  excitation  times  for  various  possible  V-T  processes  were 
obtained  from  References  15  and  16  and  plotted  in  Figure  20  as  a 
function  of  postshock  temperature.  These  theoretical  rates  are  really  ? 
smooth  curves  through  data  measured  in  shock  tubes  with  an  interferometric  J 
technique.  ' 

i 

Pure  NS  has  the  slowest  rate,  as  discussed  above,  and  it  is  obvious  > 
from  this  figure  why  it  could  not  be  examined  in  the  present  experiment.  } 

Also,  the  plot  shows  that  the  enhanced  rate  of  N!2  excitation  in  air  is 
still  quite  slow  except  at  higher  temperatures.  However,  the  O2  in  j 

air  and  pure  O2  results  demonstrate  why  O2  excitation  could  be  examined 
over  a larger  range  of  conditions.  : 

From  the  time  history  data  in  Figure  13,  the  vibrational  excitation 
time  of  oxygen  in  air  was  determined" at  T2  " 827°K.  The  .result  was 
referred  back  to  wave  fixed  coordinates  to  obtain /rn  ..1 

\ 2 "AlTy  meas  . i 

339  ysec-atm  t 11  percent.  A theoretic  value  at  tnis  temperature  was 
computed  to  be  250  nsec  using  Eqs.  (18)- (201.  The  Raman  result  is 
plotted  in  Figure  20  for  omparison  to  the  theoretical  value  for  (b  i 

excitation  in  air.  The  22  percent  variation  in  measured  excitation 
time  is  an  estimate  based  on  the  unknown  conditions  between  data  i 

points  in  Figure  13  and  the  shot  to  shot  variation  in  the  shock  speed.  I 

These  shock  speed  va;  iations  resulted  in  postshock  temperature  and  j 

excitation  time  changes  behind  the  wave.  ; 

The  measured  excitation  time  was  30  percent  slower  than  expected. 

This  is  not  an  unusual  variation  for  shock  tube  rate  data,  and  more  ; 

measurements  at  various  conditions  should  be  obtained  before  any  firm 
comments  on  accuracy  can  be  made.  However,  in  defense  of  the  Raman  ! 

data,  the  interferometric  measurements  for  in  air  are  restricted  to  ; 
temperatures  above  1500°K.  The  present  result  is  at  a temperature 
about  half  that  value  so  the  disagreement  may  be  real.  Finally,  no 
special  care  was  taken  to  assure  that  the  air  used  in  this  experiment 
was  really  dry.  Bottled  dry  air  was  used,  but  the  moisture  content  j 

was  never  checked.  Blackman's  data  showed  that  the  nrcscncc  of  water  j 

vapor  did  result  in  a slight  increase  in  the  excitation  time  as  was  • 

measured  in  this  experiment.  1 
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J 


j To  improve  the  present  excitation  time  measurement,  the  ruby  laser 

i might  possibly  be  used  in  its  "normal"  mode  instead  of  the  "Q-switched" 

| condition.  The  former  provides  a longer  pulse  (i.e.,  200-300  nsec 

j half  width).  Of  course,  the  peak  power  falls  off  accordingly  and  one 

would  pay  the  penalty  of  a loss  in  Raman  scattering  intensity  within  a 
! given  time  period.  However,  total  laser  outputs  of  10-20  joules  are 

j possible  instead  of  1-2  joules  using'  the  "Q-switch"  mode.  This 

I increased  energy  implies  that  the  Raman  scattering  output  voltage  on  , 

| the  oscilloscope  need  not  decrease  in  proportion  to  the  width  of  the 

| laser  pulse  especially  if  a suitable  RC  time  constant  was  used.  Some 
further  experiments  would  prove  useful  in  determining  if  the  "normal" 
j mode  could  be  applied  successfully.  The  biggest  advantage  to  "normal"  | 

| mode  operation  would  be  that  all  the  necessary  data  to  make  an  excitation  ; 

| time  measurement  would  be  available  from  a single  laser-shock  tube 
\ firing.  This  can  be  contrasted  to  the  17  laser  and  shock  tube  shots 
| necessary  to  obtain  a single  excitation  time  using  the  data  in 
Figure  13. 

The  temperature  and  density  ratio  measurements  obtained  from  the 
Raman  technique  for  02  and  N2  in  air  at  different  shock  speeds  are 
shown  in  Figures  21  and  22,  respectively.  The  range  of  postshock  1 

conditions  represented  by  the  data  in  tne  figures  resulted  in  a 
temperature  variation  from  7S0  to  1400*K  (i.e.,  Tj  = 300°K).  Absolute 
density  variation  cannot  be  Gxtvacted  from  Figure  22  directly,  since 
variations  in  p^  and  shock  speed  occurred.  I he  actual  postshock 
| densities  ranged  from  .35-  75  times  the  density  of  air  at  ambient 
j temperature  and  pressure  conditions.  , 

' 

Note  that  both  curve  fit  and  2-color  results,  as  described  in  the 
previous  section,  are  presented  in  Figure-  21  and  22.  The  vertical 
bars  in  the  figures  represent  the  maximum  and  minimum  ratio  for  each 
data  point  and  they  indicate  the  extremes  that  could  be  computed  due 
to  the  Raman  intensity  data  scattered  one  standard  deviation  about  the 
mean  intensity  values.  The  results  show  that  the  two  methods  of 
computing  temperature  ratios  are  in  good  agreement  with  each  other  and 
with  the  shock  tube  theory  also  plotted  in  the  figures.  The  standard 
error  is  within  12  percent  for  temperature  and  20  percent  for 
density  measurements.  This  disagreement  between  Raman  data  and  shock 
tube  theory  is  random  and  is  mostly  attributable  to  the  uncertainty 
in  the  measured  intensities.  The  uncertainty  is  mainly  due  to  the 
statistical  nature  of  the  photo-electric  effect  in  the  photomultiplier 
tube  rather  than  shot  to  shot  variation  in  shock  speed.  A greater 
quantity  of  intensity  data  would  have  improved  the  statistical  mean 
and  thereby  increased  the  accuracy  and  reduced  the  standard  deviation 
scatter.  Another  factor  that  increased  the  scatter  at  the  lower 
temperatures  was  the  decreased  sensitivity  of  the  spectra  to  a temper- 
ature change.  This  is  shown  clearly  when  one  examines  the  slope  of 
the  Tp(8260) /Ig(8285)  ratio  in  Figure  19  at  low  temperatures. 

Besides  these  factors,  several  other  experimental  effects,  which 
were  briefly  mentioned  in  the  text,  should  be  reemphasized  here  in 
the  discussion  of  the  disagreement  between  experiment  and  theory. 
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An  accurate  knowledge  of  the  actual  spectrometer  wavelength  position 
during  a specific  Raman  intensity  measurement  was  important.  The 
temperature  and  density  data  reduction  procedure  depended  upon  accurate 
A^  values.  Since  the  spectrometer  wavelength  was  changed  frequently 
between  Ag~,„  ancl  *i»  *?’  ^S*  etc*  during  a given  data  gathering  sequence 
errors  couia  have  resulted  rrom  improper  knowledge  of  the  true 
wavelength.  Of  course,  the  instrument's  indicated  wavelength  was  always 
set  by  turning  the  adjustment  knob  in  the  same  direction  to  avoid 
gear  backlash  problems  as  one  approached  Aj.  However,  hindsight 
indicates  that  a better  procedure  might0have  been  to  move  the 
wavelength  setting  to  a value  about  100A  away  from  the  desired  wave- 
length eacli  time  it  was  changed.  This  procedure  may  have  produced 
more  consistently  accurate  A values  and  Aa  settings.  The  wavelength 
inaccuracy  problem  may  be  an  explanation  for  the  larger  density 
measurement  error  compared  to  temperature  error.  In  the  density  ratio 
data  reduction  process,  not  only  was  the  basic  intensity  data  itself 
affected  by  wavelength  errors,  but  the  temperature  correction  factor 
was  also  highly  wavelength  dependent  (see  Figures  17  and  19).  Further, 
since  an  accurate  density  is  governed  by  a knowledge  of  the  correct 
temperature,  temperature  errors  compounded  the  density  measurement 
problem. 

Another  problem  was  the  change  in  sensitivity  of  the  electro- 
optical  system  with  time  due  to  instrumentation  drift  and/or  optical 
misalignment.  The  magnitude  of  this  effect  could  be  observed,  somewhat, 
by  looking  at  the  change  in  the  standard  intensity  value  with  time. 
However,  it  was  difficult  to  adjust  the  data  properly  to  account  for 
sensitivity  changes  when  purely  statistical  changes  may  have  masked 
them.  By  obtaining  all  the  data  as  quickly  as  possible,  and  allowing 
the  instrumentation  a sufficient  warm-up  time,  the  sensitivity  drift 
effect  was  minimized  but  -not  eliminated. 

The  above  factors  probably  were  the  largest  error  sources.  They 
were  much  larger  than  the  variation  in  shock  speed  from  shot  to  shot 
which  was  no  greater  than  2 percent.  Hence,  theoretical  postshock 
temperature  and  density  ratios  based  on  such  shock  speed  variation 
were  also  less  than  2 percent. 
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VII.  CONCLUSIONS 


A comparison  between  the  present  Raman  excitation  measurement 
and  data  obtained  from  interferometric  techniques  demonstrates  that 
LRS  can  be  a valuable  tool  for  determining  shock  tube  excitation  rates. 
The  inherent  ability  of  the  technique  to  measure  properties  of  individual 
gases  in  a general  gas  mixture  gives  it  an  added  advantage  over  other 
shock  tube  diagnostic  techniques,  especially  in  the  study  of  homonuclear 
molecules.  Also,  various  molecular  modes  of  vibrational  motion  in 
a poly-atomic  molecule  are  Raman  active  and  not  infrared  active. 

Perhaps  some  cases  even  exist  where  certain  molecular  modes  could 
only  be  examined  using  the  Raman  scattering  method.  In  addition  to 
V-T  excitation  times,  chemical  reaction  times,  shock  structure,  and 
other  shock  related  phenomena  could  be  examined. 

The  good  agreement  between  Raman  temperature  and  density  measure- 
ments and  known  values  obtained  from  shock  tube  theory  for 
a range  of  postshock  conditions  indicates  several  things.  The  LRS 
technique  can  independently  determine  the  density  and  temperature  of 
both  O2  and  N2  in  air  at  elevated  temperature  conditions  to  within  an 
acceptable  accuracy.  Further,  this  high  temperature  calibration 
experiment  in  conjunction  with  low  temperature  studies  in  air  and 
other  gases  implies  that  utilization  of  the  Raman  scattering  method 
should  greatly  enhance  the  amount  of  information  which  can  be  obtained 
in  general  3-dimensional  flow  fields.  The  method  could  be  applied 
to  analysis  of  the  flow  around  re-entry  vehicles,  in  rocket 
exhaust,  and  engine  plumes.  In  addition,  LRS  could  be  used  as  a 
diagnostic  technique  in  laser  cavities,  combustors,  MHD  channels 
and  wherever  else  one  desires  detailed  spatial  information  about 
the  thermochemical  state  of  a flowing  gas. 
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FIG,  3 UNRESOLVED  STOKES  ROTATION-VIBRATION  SPECTRA  OF  OXYGEN  AT  T=TR=300°K,900°  K, 
AND  1500°K  COMPUTED  FOR  THE  BROAD  EXPERIMENTAL  APPARATUS  FUNCTION. 
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FIG.  4 UNRESOLVED  STOKES  ROTATION  VIBRATION  SPECTRA  OF  OXYGEN  FOR  NONEQUILIBRIUM 
CONDITIONS  BEHIND  AN  INCIDENT  SHOCK  WAVE  FOR  TV<TR=1100*K. 
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RAMAN  SETUP  FOR  SHOCK  TUBE  EXPERIMENT  EMPHASIZING 
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FIG.  8 SCHEMATIC  OF  ELECTRO-OPTICAL  COMPONENTS  OF  DATA  GATHERING  SYSTEM 
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TIME:  50  NSEC/CM 


FIG.  9 OSCILLISCOPE  TRACE  OF  RAMAN  AND  LASER 

OUTPUT  <NZ  IN  AIR  - PA)R  = 200  TORR,  XA  = 8283A) 





FIG.  10  TIMING  CIRCUITRY  FOR  RAMAN  SHOCK-TUBE  EXPERIMENT 


NSWC/WOL/TR  75-53 


FIG.  12  DENSITY  RATIO  AS  A FUNCTION  OF  TIME  BEHIND  THE  WAVE  FRONT. 


FIG.  13  TIME  HISTORY  OF  RAMAN  SCATTERING  FROM  OXYGEN  BEHIND  A SHOCK  WAVE  IN  AIR. 
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FIG.  15  REDUCED  DATA  FOR  SCAN  METHOD  OF  DETERMINING  TEMPERATURE  AND 
DENSITY  USING  LRS. 
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Equipment 
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Naval  Intelligence 
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Intelligence  support 

Center 

4301  Suitland  Road 
Washington,  D.  C.  20390 

Department  of  Aeronautics 
DFAN 

USAF  Academy 
Colorado  80840 

COL  D.  H.  Daley,  Prof.  & 

Head 

Department  of  the  Army 
Office  of  the  Chief  of 
Research  and 
Development 
ABMDA,  The  Pentagon 
Washington,  D.  C.  20350 

Commanding  Officer 
Picatinny  Arsenal 
Dover,  New  Jersey  07801 
Mr.  A.  A.  Loeb 
SMUPA-VC-3 

Officer- in-Charge 
Naval  Weapons  Center 
Corona  Annex 

Corona,  California  91720 


3 


Commander 

y^x^cense 

New  MeK?ooe,87n59UerqUe 

Commandina  Officer 
Naval  school 

Clnei-Engineer  Cores 
Officers 

Port  Hueneme 

California  93043 

President 

Naval  War  College 

i ewport,  Rhode  island  028*0 

Superintendent 
Naval  Postgraduate 
School 
Monterey 

California  93940 
L“ing  0fficer 

Nava!  Damage  Control 

p.:faJn,inS,  Center 
Philadelphia 
Pennsylvania  19H? 

ABC  Defense  Course 


AERODYNAMICS  DEPART.MENT 
DISTRIBUTION 
(A2) 


Copies 


Copies 


Aerospace  Engineering 
Program 
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University  of 
Alabama  35486 
Prof,  W.  K.  Rey,  Chem. 

AME  Department 
University  of  Arizona 
Tucson,  Arizona  85721 
Dr,  L.  B.  Scott 

Polytechnic  Institute  of 
Brooklyn 

Graduate  Center  Library 
Route  110,  Farmingdale 
Long  Island 
New  York  11735 
Dr.  J.  Polczynski 

Polytechnic  Institute  of 
Brooklyn 
Spicer  Library 
333  Jay  Street 
Brooklyn,  New  York  11201 
Reference  Department 

Brown  University 
Division  of  Engineering 
Providence 
Rhode  Island  02912 
Dr.  M,  sibulkin 
Library 


California  Institute  of 
Technology 
Pasadena 

California  91109 

Graduate  Aeronautical 
Laboratories  Aero. 
Librarian 
Dr . H . Liepmann 
Karman  Lab-301 
Prof.  L.  Lees,  Firestone 
Flight  Science  Lab. 

Dr.  D.  Coles, 

306  Karman  Lab, 

Dr.  A.  Roshko 

University  of  California 
Berkeley 

California  94720 

Dr.  M.  Holt,  Div.  of 
Aeronautical  Sciences 
Research  Coordinator 
College  of  Engineering 

GAS DYNAMICS 

University  of  California 
Richmond  Field  Station 
1301  South  46th  Street 
Richmond 

California  94804 
A,  K.  Oppenheim 

Department  of  Aerospace 
Engineering 
University  of  Southern 
California 
University  Park 
Los  Angeles 
California  90007 
Dr.  John  Laufer 
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University  of  California  - 
San  Diego 

Department  of  Aerospace 
and  Mechanical 
Engineering  Sciences 
LaJolla 

California  92037 
Dr,  P.  A.  Libby 

Case  Western  Reserve 
University 
3306  Clarendon  Road 
Cleveland,  Ohio  44118 
G,  Kuerti 

The  Catholic  University 
of  America 

Washington,  D.  C.  20017 
Dr,  C,  C.  Chang 
Dr,  Paul  K.  Chang 

Mechanical  Engr.  Dept. 
Dr.  M.  J.  Casarella 

Mechanical  Engr.  Dept. 

University  of  Cincinnati 
Cincinnati,  Ohio  45221 
Department  of  Aerospace 
Engineering 
Dr,  Arnold  Polak 

Department  of  Aerospace 
Engineering  Sciences 
University  of  Colorado 
Boulder,  Colorado  80302 

Cornell  University 
Graduate  School  of  Aero. 
Engineering 

Ithaca,  New  York  14850 
Prof,  W.  R,  Sears 
Dr,  S.  F.  Shen 
Prof.  F.  K.  Moore,  Head 
Thermal  Engineering 
Dept.,  208  Upson  Hall 

University  of  Delaware 
Mechanical  and  Aeronautical 
Engineering  Dept, 

Newark , Delaware  19711 
Dr,  James  E,  Danberg 


Georgia  Institute  of 
Technology 

225  North  Avenue,  N.W. 

Atlanta,  Georgia  30332 
Dr,  Arnold  L,  Ducoffe 
Aerospace  Engineering 

Dept.  , 

Technical  Reports  Collection 
Gordon  Mckay  Library 

Harvard  University  | 

Div.  of  Eng'g.  and  Applied  ; 

Physics 

Fierce  Hall,  Oxford  Street 
Cambridge 

Massachusetts  02138 

Illinois  Institute  of 
Technology 

3300  South  Federal  1 

Chicago,  Illinois  60616 
Dr.  M.  V.  Morkovin 
Prof.  A.  A.  Fejer, 

M.A.E.  Dept. 

University  of  Illinois  , 

101  Transportation  Building 
Urbana,  Illinois  61801 
Aeronautical  and  Astro- 
nautical  Engineering 
Dept. 

Iowa  State  University 
Ames,  Iowa  50010 

Aerospace  Engineering 
Dept. 

The  Johns  Hopkins  University 
Baltimore,  Maryland  21218 
Prof.  S.  Corrsin 

University  of  Kentucky 
Wenner-Gren  Aero,  Lab. 

Lexington,  Kentucky  40506 
C.  F.  Knapp 


2 


Conies 


Copies 


Department  of  Aero, 
Engineering  ME  106 
Louisiana  State 
Unive3  sity 
Baton  Rouge 
Louisiana  70803 
Dr,  R,  H.  Miller 

University  of  Maryland 
College  Park 
Maryland  20740 

Prof,  A.  Wiley  Sherwood 
Department  of  Aerospace 
Engineering 

Prof.  Charles  A.  Shreeve 
Department  of  Mechanical 
Engineering 

Dr,  S.  I.  Pai,  Institute 
of  Fluid  Dynamics 
and  Applied 
Mathematics 
Dr,  Redfield  W.  Allen 
Department  of 
Mechanical  Engineering 
Dr,  W.  L.  Melnik 
Department  of 
Aerospace  Engr. 

Dr,  John  D.  Anderson,  Jr. 
Department  of 
Aerospace  Engineering 


Michigan  State 

fni\?ersity  Library 
East  Lansing,  Michigan 
Documents  Department 


48823 


Massachusetts  Institute  of 
Technology 
Cambridge 

Massachusetts  02139 
Mr,  J.  R.  Martuccelli 
Rm  33-211 
Prof.  M.  Finston 
Prof.  J,  Baron,  Dept,  of 
Aero,  and  Astro., 

Rm,  37-461 

Prof,  A.  H.  Shapiro,  Head 
Mech.  Engr,  Dept. 

Aero.  Engineering  Library 
Prof,  Ronald  F,  Probestein 
Dr.  E.  E,  Covert 

Aerophvsics  Laboratory 


University  of  Michigan 
Ann  Arbor,  Michigan  48104 
Dr,  A.  Kuethe,  Dept,  of 
Aero.  Engineering 
Dr.  M.  Sichel,  Dept. 

of  Aero.  Engineering 
Engineering  Library 
Aerospace  Engineering 
Library  Aquisitions 
Section  Serials 
Division 

Serials  and  Documents 
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General  Library 
University  of  Michigan 
Ann  Arbor,  Michigan  48104 

Mississippi  State 
University 

Department  of  Aerophysics 
and  Aerospace  Engineering 
Drawer  AP 

Mississippi  39762 

Mr,  Charles  L.  Cliett 

U,  S.  Naval  Academy 
Annapolis,  Maryland  21402 
Engineering  Department, 
Aerospace  Division 

Library,  Code  2124 
U.  S.  Naval  Postgraduate 
School 
Monterey 

California  93940 
Technical  Reports 
Section 

New  York  University 
University  Heights 
New  York,  New  York  10453 
Dr,  Antonio  Ferri,  Director 
of  Guggenheim  Aerospace 
Laboratories 
Prof,  V,  Zakkay 
Engineering  and  Science 
Library 


Copies 


Copies 


North  Carolina  State 
College 
Raleigh 

North  Carolina  276C7 
Dr,  R,  W.  Truitt,  Head, 
Mech , and  Aero . 
Engineering 

Dr,  H.  A.  Hassan,  Dept, 
of  Hech.  and  Aero, 
Engineering 

D.  H,  Hill  Library 
North  Carolina  State 
University 
P,  0,  Box  5007 
Raleigh 

North  Carolina  27607 

University  of  North 
Carolina 
Chapel  Hill 
North  Carolina  27514 
Department  of  Aero. 
Engineering 

Library,  Documents  Section 
AFROTC  Det  590 

Northwestern  University 
Technological  Institute 
Evanston,  Illinois  60201 
Department  of  Mechanical 
Engineering 
Library 

Dr.  John  D.  Nicolaides 
P.  0.  Box  556 

Notre  Dame,  Indiana  46556 

Ohio  State  University 
2036  Neil  Avenue 
Columbus,  Ohio  43210 
Aero.  Civil  Library 
Prof.  J.  D.  Lee 

Aeronautical  Research 
Lab . 

Prof,  G,  L,  Von  Eschen 
Dept,  of  Aero-Astro 
Engr, 


Ohio  State  University 
Libraries 

Documents  Division 
1858  Neil  Avenue 
Columbus,  Ohio  43210 

Department  of  Aerospace 
Engineering 

Room  233  Hammond  Building 
The  Pennsylvania  State 
University 
University  Park 
Pennsylvania  16802 

Pennsylvania  State 
University  Library 
Documents  Section 
University  Park 
Pennsylvania  16802 

Bevier  Engineering 
Library 

126  Benedum  Hall 
University  of 
Pittsburgh 
Pittsburgh 
Pennsylvania  15261 
Librarian 

Princeton  University 
James  Forrestal  Research 
Center 

Gas  Dynamics  Laboratory 
Princeton,  New  Jersey  08540 
Prof.  S.  Bogdonoff 
Mr.  I.  E,  Vas 

Purdue  University 
School  of  Aeronautical 
and  Engineering  Sciences 
Lafayette,  Indiana  47907 
Library 

Dr.  P,  S,  Lykoudis,  Dept, 
of  Aero,  Engineering 

Rensselaer  Polytechnic 
Institute 

Troy,  New  York  12181 
Dept,  of  Aeronautical 
Engineering  and 
Astronautics 
Dr,  Robert  E,  Duffy 
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Rutgers  - The  State 
University 

University  Heights  Campus 
New  Brunswick 
New  Jersey  08903 

Dr.  R.  H.  Page,  Dept, 
of  Mech.  and  Aero, 
Engineering  2 

Stanford  University 
Stanford 

California  94305 
Librarian,  Dept,  of 
Aeronautics  and 
Astronautics 

Stevens  Institute  of 
Technology 

Hoboken,  New  Jersey  07030 
Mechanical  Engineering 
Department 
Library 

The  University  of  Texas 
at  Austin 
Applied  Research 
Laboratories 
P.  0.  Box  8029 
Austin,  Texas  78712 
Director 

University  of  Toledo 
Department  of  Aero. 

Engineering 
Research  Foundation 
Toledo,  Ohio  43606 

Documents  Department 
Virginia  Polytechnic 
Institute 
Blacksburg 
Virginia  24061 
Carol  M.  Newman 
Library 


University  of  Virginia 
Charlottesville 
Virginia  22901 
Alderman  Library 
Science  Reference 
Division 

Dr.  G.  Matthews,  Dept, 
of  Aerospace 
Engineering 

University  of  Washington 
Seattle,  Washington  98105 
Engineering  Library 
Dept,  of  Aeronautics 
and  Astronautics 
Prof,  k.  E.  Street,  Dept. 

of  Aero,  and  Astro. 

Prof.  A.  Hertzberg,  Aero, 
and  Astro.  Guggenheim 
Hall 

West  Virginia  University 
Morgantown 
West  Virginia  26506 
Library 

Federal  Reports  Center 
University  of  Wisconsin 
462  Mechanical  Engineering 
Building 

1513  University  Avenue 
Madison,  Wisconsin  53706 
S,  Reilly 
J.  S,  Murphy 

Los  Alamos  Scientific 
Laboratory 
P,  0,  Box  1663 

Los  Alamos,  New  Mexico  87544 
Report  Library 
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University  of  Maryland 
Baltimore  County  (UMBC) 

5401  Wilkens  Avenue 
Baltimore 
Maryland  21228 
Dr.  R.  C.  Roberts 
Mathematics  Dept. 

Oklahoma  State  University 
Office  of  Engineering 
Research 
Stillwater 
Oklahoma  74074 

Dr.  V.  S.  Haneman,  Jr. 

Institute  for  Defense 
Analyses 

400  Army-Navy  Drive 
Arlington 
Virginia  22202 

Classified  Library 

Kaman  Nuclear 
1700  Garden  of  the 
Gods  Road 
Colorado  Springs 
Colorado  80907 

Kaman  Science  Corporation 
Avidyne  Division 
83  Second  Avenue 
Burlington 
Massachusetts  01803 
Dr.  J.  R.  Ruetenik 

Rockwell  International 
B-l  Division 
International  Airport 
Technical  Information  Center 
(BAGS ) 

Los  Angeles,  California  90009 
LAD  Library,  Dept.  299 


North  American  Rockwell 
Corporation 

Engineering  Data  Services 
4300  E.  Fifth  Avenue 
Columbus,  Ohio  43216 

M.I.T.  Lincoln  Laboratory 

P.  0.  Box  73 

Lexington 

Massachusetts  02173 
Library  A-082 
Dr.  S.  Edelberg 
Dr.  R.  H.  Kingston 
Mr.  J.  Freedman 
Dr.  G.  P.  Dinneen 
Dr.  R.  H.  Rediker 

The  RAND  Corporation 
1700  Main  Street 
Santa  Monica 
California  90406 
Library  - D 

The  Boeing  Company 
P.  0.  Box  3999 
Seattle 

Washington  98124 
Aerospace  Library 

8K-38  J.  M.  MacDonald 

United  Aircraft  Corporation 
Research  Laboratories 
East  Hartford 
Connecticut  06108 
Dr.  William  M.  Foley 

United  Aircraft  Corporation 
400  Main  Street 
East  Hartford 
Connecticut  06108 
Library 

Hughes  Aircraft  Company 
Continela  and  Teale  Streets 
Culver  City 
California  90230 

Co.  Tech.  Doc.  Ctr. , 
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Lockheed  Missiles  and 
Space  Company 
32 si  Hanover 
Palo  Alto 
California  9^304 
H.  S.  Bettencourt 
Dent.  81-11 
Bldg.  154 
L.  R.  Lunsford 
Orgn.  52-24 
Bldg.  201 

Lockheed  Missiles  and 
Space  Company 
3251  Hanover  Street 
Palo  Alro 
California  94304 

Technical  Information 
Center 

Lockheed-California 

Company 

Burbank 

California  91503 
Central  Library 
Dept.  84-40 
Bldg.  170,  PLT.  B-l 

Vice  President  and 
Chief  Scientist 
Dept.  03-01 
Lockheed  Aircraft 
Corporation 
P.  0.  Box  551 
Burbank 

California  91520 
Dr.  Ronald  Smelt 

Martin  Marietta  Corporation 
P.  0.  Box  988 
Baltimore 
Maryland  21203 

Science-Technology 

Library  {Mail  No,  398) 


Martin  Company 
3211  Trade  Winds  Trail 
Orlando,  Florida  32805 
Mr.  H.  J.  Diebolt 

General  Dynamics 
P.  0.  Box  748 
Fort  Worth,  Texas  76101 
Research  Library 
2246 

George  Kaler,  Mail 
Zone  2880 

Cornell  Aeronautical 
Laboratory,  Inc. 

4455  Genesee  Street 
Buffalo,  New  York  14221 
Library 

Air  University  Library 
(SE)  63-578 

Maxwell  Air  Force  Base 
Alabama  36112 

McDonnell  Douglas 
Corporation 
Research  & Engineering 
Library 

Dept.  209,  Bldg.  33 
P.  0.  Box  516 
3t.  Louis 
Missouri  63166 

McDonnell  Douglas 

Astronautics  Company 
Space  Systems  Center 
5301  Bolsa  Avenue 
Huntington  Beach 
California  90405 
A2-260  Library 
Dr.  J,  S.  Murphy,  A-830 
Mr.  W,  H.  Branch,  Director 
A3-339  Library 
P.  L.  Klenatt" A3-833 
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Fairchild  Hiller 
Republic  Aviation 
Division 
Farmingdale 
New  York  11735 

Engineering  Library 

General  Applied  Science 
Laboratories,  Inc. 

Merrick  and  Stewart 
Avenues 

Westburv,  Long  Island 
New  York  11590 
Dr.  F.  Lane 

The  Whitney  Library 
General  Electric  Research 
and  Development  Center 
The  Knolls,  '1-1 
P.  0.  Box  8 
Schenectady 
New  York  12301 

M.  F.  Orr,  Manager 

General  Electric  Company 
Missile  and  Space 
Division 
P.  0.  Box  8555 
Philadelphia 
Pennsylvania  19101 
MSD  Library 

Larry  Chasen,  Mgr. 

Dr.  J.  D.  Stewart,  Mgr. 

Research  and  Engineering 

General  Electric  Company 
AEG  Technical  Information 
Center  N-32 

Cincinnati,  Ohio  45215 


Copies 

General  Electric  Company 
Missile  and  Space 
Division 
P.  0.  Box  8555 
Philadelphia 
Pennsylvania  19101 
Dr.  S.  M.  Scala 
Dr . H . Lew 
Mr.  J.  W.  Faust 
Mr.  W.  Daskin 
S.  B.  Kottock 
S . Kahn 
J.  B.  Arnaiz 
L.  A.  Marshall 
A.  Martellucci 

AVCO-Everett  Research 
Laboratory 

2385  Revere  Beach  ^arkway 
Everett 

Massachusetts  02149 
Library 

Dr.  George  Sutton 

LTV  Aerospace  Corporation 
Vought  Systems  Division 
P,  O.  Box  5907 
Dallas,  Texas  75222 

Mr.  F.  G.  Simpson  2 

Mail  Station  2-54142 

LTV  Aerospace  Corporation 
Missiles  and  Space  Division 
P.  0.  Box  6267~ 

Dallas,  Texas  75222 
MSD-T-Library 

Northrop  Norair 
3901  West  Broadway 
Hawthorne 
California  90250 
Tech.  Info.  3343-32 

Grumman  Aircraft  Engineering 
Corporation 
Bethpage,  Long  Island 
New  York  11714 

Mr.  R,  A.  Scheuing 
Mr,  H.  B,  Hopkins 
Mr.  H.  R.  Reed 
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Marquardt  Air'  raft 
Corporation 
16555  Saticoy  Street 
Van  Nuys 

California  91409 
Library 

ARDE  Associates 
P.  0,  Bos  286 
580  Winters  Avenue 
Paramus,  New  Jersey  07652 
Librarian 

Aeronautical  Research 
Associates  of 
Princeton 

50  Washington  Road 
Princeton , New  Jersey  08540 
Dr.  C.  duP.  Donaldson 

General  Research  Corporation 
5383  Hollister  Avenue 
P.  0.  Box  3587 
Santa  Barbara 
California  93105 
Technical  Info.  Officer 

Sandia  Corporation 
Sandia  Base 
Albuquerque 
New  Mexico  87115 
Mr.  K.  Coin,  Div.  9322 
Mrs.  B.  R.  Allen 
Div  3421 

HERCULES  INCORPORATED 
Allegany  Ballistics 
Laboratory 
P.  0.  Box  210 
Cumberland 
Maryland  21502 
Mrs.  Louise  E.  Derrick 
Librarian 

General  Electric  Company 
P.  O.  Box  2500 
Daytona  Beach 
Florida  32015 
Dave  Ho  vis , i*n.  4109 


TRW  Systems  Group 
1 Space  Park 
Redondo  Beach 
California  90278 
Technical  Library- 
Document 
Acquisitions 

Stanford  Research 
Institute 

333  Ravenswood  Avenue 
Menlo  Park 
California  94025 
Dr.  G.  Abrahamson 
Mr.  John  Malick  L3088 

Hughes  Aircraft  Company 
P.  0.  Box  3310 

Bldg.  .00/MSE-245 
Fullerton 
California  92634 

R.  H.  Sterling,  600-E201 

Westinghouse  Electric 
Corporation 

Astronuclear  Laboratory 
P.  0.  Box  10864 
Pittsburgh 
Pennsylvania  15236 
Library 

CONVAIR  Division  of  General 
Dynamics 

Library  and  Information 
Services 
P.  0,  Box  12009 
San  Diego 
California  92112 

CONVAIR  Division  of 
General  Dynamics 
P.  0.  Box  1128 
San  Diego 
California  921.1,2 
Dr.  Jan  Raat^  Aero- 
ballistics 

Dept.  Mail  Zone  583-00 
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AVCO  Missiles  Systems 
Division 

201  Lowell  Street 
Wilmington 
Massachusetts  01887 
E.  E.  H,  Schumann 
J.  Otis 

Chrysler  Corporation 
Space  Division 
P.  O.  Box  29200 
New  Orleans 
Louisiana  70129 

G.  T.  Boyd,  Dept,  2781 

General  Dynamics 
Pomona  Division 
Pomona,  California  91766 
Division  Library 
Mail  Zone  6-20 

Phil co-Ford  Corporation 
Aeroneutronic  Division 
Newport  Beach 
California  92660 
Dr.  A.  Demetriades 
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Raytheon  Company 
Missile  Systems  Division 
Hartwell  Hoad 
Bedford 

Massachusetts  01730 
Dr,  H.  A,  Mehlhom 
Optical  Systems 
Dept.  (S4-55) 

TRW  Systems  Group 
Space  Park  Drive 
Houston,  Texas  77058 
M.  W.  Sweeney,  Jr. 

Marine  Bioscience 
Laboratory 
527  Las  Alturas  toad 
Santa  Barbara 
California  93103 
Dr.  A.  C.  Charters 


University  of  Wyoming 
Laramie,  Wyoming  82070 
Library 
Prof.  Matheny 

ELMAC,  Division  Varian 
Associates 
301  Industrial  Way 
San  Carlos 
California  S4070 

J.  R.  Quinn,  Marketing 
Manager 

The  Aerospace  Corporation 
P.  Q.  Box  92957 
Los  Angeles 
California  90009 

Dr.  W.  R.  Warren,  Jr. 
130/691 

Department  of  the  Air  Force 
Headquarters  Air  Force 
Special  Communications 
Center 
San  Antonio 
Texas  78243 
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Defense  Advanced 
Research  Projects 
Agency 

1400  Wilson  Boulevard 

Arlington 

Virginia  22209 

Director,  Laser  Division 

ODDR&E 

Pentagon 

Washington,  D,  C,  20301 
Assistant  Director 
(Space  and  Advanced 
Systems) 

U.  S.  Arms  Control  and 
Disarmament  Agency 
Department  of  State 
Building,  Room  4931 
Washington,  D.  C.  20451 
Dr.  Charles  Henkin 

National  Aeronautics  and 
Space  Administration 
Code  RR,  ROB  10B 
600  Independence  Avenue,  S.W. 
Washington,  D.  C.  20546 

Division  of  Military 
Applications 

Energy  Research  and  Development 
Administration 

Washington,  D.  C.  20545 
Dr.  Lawrence  Killion 

National  Aeronautics  and 
Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio  44135 
Dr.  John  W.  Dunning,  Jr. 
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